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CHAPTER | 


INTRODUCTION 


The term FORTRAN IV (FORmula TRANslation) is used interchangeably to designate both the FORTRAN IV 
language and the FORTRAN IV translator or compiler. The FORTRAN IV language is composed of 
mathematical-form statements constructed in accordance with precisely formulated rules. FORTRAN IV 
programs consist of meaningful sequences of FORTRAN statements intended to direct the computer to 


perform the specified operations and computations. 


The FORTRAN IV compiler is itself a computer program that examines FORTRAN IV statements and tells 
the computer how to translate the statements into machine language. The compiler runs in a minimum of 
9K of core. The program written in FORTRAN IV language is called the source program. The resultant 


machine language program is called the object program. 


FORTRAN IV includes such advanced features as logical operators, type declaration statements, double 


precision and complex arithmetic, named COMMON, and DATA statements. 


FORTRAN IV language elements are discussed in Chapter 2 of this manual, followed by separate chapters 
on the five categories of FORTRAN IV statements (arithmetic, control, input/output, specification, and 
subprogram). The appendices contain a list of FORTRAN statements and summary descriptions of library 


functions and subroutines. 


LINE FORMAT 


Each line of a FORTRAN program consists of three fields: statement number field, line continuation field, 


and statement field. 


Statement Number Field 


A statement number consists of from one to five digits. Leading zeros and all blanks in 
this field are ignored. Statement numbers may be in any order and must be unique. Any statement ref- 


erenced by another statement must-have a statement number . 


Line Continuation Field 
If a FORTRAN statement is so large that it cannot conveniently 
fit into one statement field, the statement fields of up to 19 
additional lines may be used to specify the complete statement. 
Any line which is not continued or the first line of a continued 
statement must have a blank or zero in Column 6 ‘see "use of tab"). 
Continuation lines must have a character other than blank or zero 
in Column 6. 

Use of Tab 

A horizontal tab may be used in any of the following 3 ways: 


1. <TAB> "statement field" 


2. "statement number" <TAB> "Statement field" 
3. <TAB> "number" "statement field" 
"number" = digits 1 thru 9 used to denote line continuation. 


Statement Field 
Any FORTRAN statement, as described in later sections, may appear in 
the statement field (columns 7-72). Except within the literal 
descriptor of a FORMAT statement, blanks (or spaces) are ignored and 


may be used freely for appearance purposes. 


Comment Lines 
Any line which starts with the letter C in column 1 is interpreted 
as a line of comments. Comment lines are printed onto any listings 
requested but are otherwise ignored by the compiler. Columns 2-72 


may be used in any format for comment purposes. 


CHARACTER SET 


The following characters are used in the FORTRAN IV language: 


Blank 0 @ 
2 1 A 
ty 2 B 
# 3 Cc 
$ 4 D 
% 5 E 
& 6 F 
CY) 7 G 
( 8 H 
) 9 ! 

% : J 
+ K 
; < L 
- = M 
: > N 
/ ? Oo 


The following characters are not used: C\NI+ 


~N<xXs<CANzH? 


CHAPTER 2 


CONSTANTS, VARIABLES, AND EXPRESSIONS 


The rules for defining constants and variables and for forming expressions to evaluate functions are described 


in this chapter. 


CONSTANTS e 


Seven types of constants are permitted in a FORTRAN IV source program: integer or fixed point, real or 


single-precision floating point, double-precision floating point, octal, complex, logical, and literal. . 


Integer Constants 


An integer constant consists of from one to eleven decimal digits written without a decimal point. 


EXAMPLES: 3 
-528 
8085 
: rere 35 35 ; 
An integer constant must fall within the range - 2+] to 2. -1. When used for the value of a subscript 


8 


: ; 1 
or as an index in a DO statement, the value of the integer is taken as modulo 2. 


Real Constants 


Real constants are written as a string of decimal digits including a decimal point. Any number of digits, 
of which nine are significant, may be written. Real constants may be given a decimal scale factor by 
appending an E followed by an integer constant. The field following the letter E must not be blank, but 


t 
may be zero. 


EXAMPLES: 15. 
.579 
5.0E3(i.e., 5000.) 


A real constant has precision to eight digits. The magnitude must lie approximately within the range 


0.14E-38 to 1.7E38. 


1 


Double Precision Constants 


A double precision constant is specified by a string of decimal digits, including a decimal point, which 
are followed by the letter D and the decimal scale factor. The field following the letter D must not be 


blank, but may be zero. 


EXAMPLES: 24.671 325982134D0 
3.6D2 (i.e., 360.) 
3.6D-2 (i.e., .036) 


Double precision constants have precision to 16 digits. The magnitude of a double precision constant 


must lie approximately between 0.14E-38 and 1.7E38. 


Octal Constants 


A number preceded by a double quote represents an octal constant. An octal constant may appear in an 
arithmetic or logical expression or a DATA statement. Only the digits 0-7 may be used and only the first 


twelve digits are significant. 


EXAMPLES: "7777 
"-31563 


Complex Constants 


FORTRAN IV provides for direct operations on complex numbers. Complex constants are written as an 


ordered pair of real constants separated by a comma and enclosed in parentheses. 


EXAMPLES: (.70712, -.70712) 
(8.763E3, 2.297) 


The first constant of the pair represents the real part of the complex number, and the second constant 
represents the imaginary part. The real and imaginary parts may each be signed. The enclosing paren- 
theses are part of the constant and always appear, regardless of context. 


If A= a, tia, and B = b, 


become: 


+ib, (where i = V-1), FORTRAN IV arithmetic operations on complex numbers 


AtB = a,b, +i(a,tb,) 


AsB = (ab, -a,b,)+i (a,b) +0)b.) 


ve seal ar a (a,b, -a)b,) 
2, 2 2 2 
b, +b, b, +b, 


Logical Constants 


The two logical constants, .TRUE. and .FALSE., have the internal values—1 and 0, respectively. The 


enclosing periods are part of the constant and always appear. 


Logical constants may be entered in DATA or input statements as signed octal integers (—1 and 0). Logical 
quantities may be operated on in either arithmetic or logical statements. Only the sign is tested to de- 


termine the truth value of a logical variable. 


Literal Constants 


A literal constant may be in either of two forms: 


1. Astring of characters enclosed in single quotes; two adjacent single quotes within 


the constant are treated as one single quote. 


2. Astring of the form: 


nHx,)X.- Xe 


where x)X5-- +X) is the constant, and n is the number of characters following the H. 


EXAMPLES: ‘LITERAL CONSTANT! 
‘DON''T! 
5HDON'T 


VARIABLES 


A variable is a quantity whose value may change during the execution of a program. Variables are speci- 
fied by name and type. The name of a variable consists of one or more alphanumeric characters, the first 
one of which must be alphabetic. Only the first six characters are interpreted as defining the variable 
name. The type of variable (integer, real, logical, double precision, or complex) may be specified by 

a type declaration statement, described in Chapter 6. Variables of any type may be either scalar or 


array variables. 


SCALAR VARIABLES 


A scalar variable represents a single quantity. 


EXAMPLES: A 
G2 
POPULATION 


ARRAY VARIABLES 


An array variable represents a single element of a one-to-n dimensional array of quantities. The variable 
is denoted by the array name followed by a subscript list enclosed in parentheses. The subscript list is a 
sequence of integer expressions, separated by commas. The expression may be arithmetic combinations 

of integer variables or integer constants. Each expression represents a subscript, and the values of the 
expressions determine the array element referred to. For example, the row vector A. would be represented 
by the subscripted variable A(J), and the element, in the second column of the first row of the square ma- 


trix A, would be represented by A(1,2). Arrays may have any number of dimensions. 


EXAMPLES: Y(1) 
STATION (K) 
A (3* K+2, I, J-1) 


The three arrays above (Y, STATION, and A) would have to be dimensioned by a DIMENSION, COM- 
MON, or type declaration statement prior to their first appearance in an executable statement or ina 


DATA or NAMELIST statement. (Array dimensioning is discussed in chapter 6.) 


Arrays are stored in increasing storage locations with the first subscript varying most rapidly and the last 


ve 


subscript varying least rapidly. For example, the 2-dimensional array B(I,J) is stored in the following 
order: B(1,1), B(2,1),..-, B(I,1),B (1,2),B (2,2),...,B (I,2),...,B (I,J). 


EXPRESSIONS 


Expressions may be either numeric or logical. To evaluate an expression, the object program performs 


the calculations specified by the quantities and operators within the expression. 


Numeric Expressions 


A numeric expression is a sequence of constants, variables, and function references separated by numeric 


operators and parentheses in accordance with mathematical convention and the rules given below. 


The numeric operators are+, —, *, /, **, denoting, respectively, addition, subtraction, multiplication, 


division, and exponentiation. 


In addition to the basic numeric operators, function references are also provided to facilitate the evaluation 
of functions such as sine, cosine, and square root. A function is a subprogram which acts upon one or more 
quantities, called arguments, to produce a single quantity called the function value. Function references 
are denoted by the identifier, which names the function (such as SIN, COS, etc.), followed by an argu- 


ment list enclosed in parentheses: 
identifier(argument, argument, ..., argument) 
At least one argument must be present. An argument may be an expression, an array identifier, a subpro- 


gram identifier, or an alphanumeric string. 


Function type is given by the type of the identifier which names the function. The type of the function 


is independent of the types of its arguments. 
A numeric expression may consist of a single element (constant, variable, or function reference): 


2.71828 
Z(N) 
TAN(THETA) 


Compound numeric expressions may be formed by using numeric operators to combine basic elements: 


X+3. 
TOTAL/A 
TAN(PI*M) 


* 


Compound numeric expressions must be constructed according to the following rules: 


1. With respect to the numeric operators +, -, *, /, ony type of quantity (logical, 
octal, integer, real, double precision, complex or litercl) may be combined with any 
other, with one exception: a complex quantity cannot be combined with a double 


precision quantity. 


The resultant type of the combination of any two types may be found in Table 1. The 


conversions between data types will occur as follows: 


(a) A literal constant will be combined with any other integer constant as an 
integer and with a real or double word as a real or double word quantity. 


(Double word refers to both double precision and complex.) 


(b) An integer quantity (constant, variable, or function reference) combined 
with a real or double word quantity results in an expression of the type real 

or double word respectively; e.g., an integer variable plus a complex variable 
will result in a complex subexpression. The integer is converted to floating 
point and then added to the real part of the complex number. The imaginary 


part is unchanged. 


(c) Areal quantity (constant, variable, or function reference combined with 
a double word quantity results in an expression that is of the same type as the 


double word quantity. 


(d) A logical or octal quantity is combined with an integer, real, or double 


word quantity as if it were an integer quantity in the integer case, or a real 


quantity in the real or double word case (i.e., no conversion takes place). 


2. Any numeric expression may be enclosed in parentheses and considered to be a basic 


element. 
(X+Y)/2 
(ZETA) 
(COS(SIN(PleM} X)) 


TABLE 1 TYPES OF RESULTANT SUBEXPRESSIONS 


Type of Quantity 


Logical 
Double ! 
eee Complex nes 
Complex Double Real 
Precision 
® 
Complex Double Integer 
Precision 


Complex Complex Complex Complex Not Complex 
Allowed 


Double Double Not Double Double 
Precision Precision Allowed Precision Precision 


Double 
Precision 


Type of Quantity 


Logical, Complex Double Logical, 
Octal, or Precision Octal, or 
Literal Literal 


3. Numeric expressions which are preceded by a + or — sign are also numeric expressions: 


+X 
-(ALPHA*BETA) 
-SQRT(-GAMMA) 


4. If the precedence of numeric operations is not given explicitly by parentheses, it 


is understood to be the following (in order of decreasing precedence): 


Operator . 
a numeric exponentiation 

*and/ numeric multiplication and division ‘ 
+and- numeric addition and subtraction 


In the case of operations of equal hierarchy, the calculation is performed from left 


to right. This is also true for exponentiation. 


t See also page 14 


5. No two numeric operators may appear in sequence. For instance: 
X*«-Y 
is improper. Use of parentheses yields the correct form: 
X«(-Y) 
By use of the foregoing rules, all permissible numeric expressions may be formed. As 


an example of a typical numeric expression using numeric operators and a function 


reference, the expression for the largest root of the general quadratic equation: 


i be Sidbe 


2a 
would be coded as: 


(—B+SQRT (B¥*#2—4.*A*C))/(2.*A) 


Logical Expressions 


A logical expression consists of logical constants, logical variables, logical function references, and 


arithmetic expressions, separated by logical operators or relational operators. Logical expressions are 


provided in FORTRAN IV to permit the implementation of various forms of symbolic logic. Logical con- 
stants are defined by arithmetic statements, which are described in Chapter 3. Logical variables and 
functions are defined by the LOGICAL statement, described in Chapter 6. Binary variables may be rep- 


resented by the logical constants .TRUE. and .FALSE., which must always be written with enclosing periods. 


Logical Operators 


The logical operators, which include the enclosing periods and their definitions, are as follows, where P 


and Q are logical expressions: 


-.NOT.P Has the value .TRUE. only if P is .FALSE., and has the value 
.FALSE. only if P is . TRUE. 


P.AND.Q Has the value . TRUE. only if P and Q are both .TRUE., and 
has the value .FALSE. if either P or Q is . FALSE. 


P.OR.Q (Inclusive OR) Has the value .TRUE. if either P or Q is .TRUE., 
and has the value .FALSE. only if both P and Q are .FALSE. 


P.XOR.Q (Exclusive OR) Has the value . TRUE. if either P or Q but not ? | 
both are .TRUE., and has the value .FALSE. otherwise. 


P.EQV.Q (Equivalence) Has the value . TRUE. if P and Q are both 
. TRUE. or both .FALSE., and has the value .FALSE. otherwise. 


Relational Operators 


The relational operators are as follows: 


Operator Relation 
.GT. greater than 
.GE. greater than or equal to 
kT less than 
“LE less than or equal to 
.EQ. equal to 
-NE. not equal to 


The enclosing periods are part of the operator and must be present. 


Mixed expressions involving integer, real, and double precision types may be combined with relationals. 


The value (.TRUE. or . FALSE.) of such relations will be calculated by subtraction; i.e., 
expression, "relation" expression, 
will be calculated as though: 
expression, - expression. "relation" zero 
had been written. 
The relational operators .EQ. and .NE. may also be used with COMPLEX expressions. (Double word 
quantities are equal if the corresponding parts are equal . ) 


A logical expression may consist of a single element (constant, variable, function reference, or relation): 


- TRUE. 
X.GE.3.14159 


Single elements may be combined through use of logical operctors to form compound logical expressions, 
such as: 


TVAL.AND.INDEX 
BOOL(M).OR.K.EQ.LIMIT 


Any logical expression may be enclosed in parentheses and regarded as an element: 


(T.XOR.S).AND.(R.EQV.Q) 
PARITY ((2.GT.Y.OR.X.GE.Y). AND. NEVER) 


Any logical expression may be preceded by the unary operator. NOT. as in: 


-NOT.T 
»~NOT.X+7.GR.Y+Z 
BOOL(K). AND. .NOT.(TVAL .OR.R) 


No two logical operators may appear in sequence, except in the case where .NOT. appears as the second 


of two logical operators, as in the example above. 


Two decimal points may appear in sequence, as in the example above, or when one belongs to an operator 


and the other to a constant. 


When the precedence of operators is not given explicitly by parentheses, it is understood to be as follows 


(in order of decreasing precedence): 


kk 


«,/ 

tiie 
.GT.,.GE.,.LT.,.LE.,.EQ.,.NE. 
.NOT. 

.AND. 

OR. 

-EQV., .XOR. 


For example, the logical expression 
-NOT. ZETA**2+Y*MASS.GT.K—-2.OR.PARITY .AND.X.EQ.Y 


is interpreted as 


(.NOT. (((ZETA#*2)}+ (Y*MASS)). GT. (K—2))). OR. (PARITY . AND. (X.EQ.Y)) 


ek: ae 


CHAPTER 3 


THE ARITHMETIC STATEMENT 


One of the key features of FORTRAN IV is the ease with which arithmetical computations can be coded. 
Computations to be performed by FORTRAN IV are ind’ :ated by arithmetic statements, which have the 


general form: 
A=B 


where A is a variable, B is an expression, and = is a replacement operator. The arithmetic statement 
causes the FORTRAN IV object program to evaluate the expression B and assign the resultant value to the 


variable A. Note that the = sign signifies replacement, not equality. Thus, expressions of the form: 


A=A+B and 
A=AxB 


are quite meaningful and indicate that the value of the variable A is to be changed. 


EXAMPLES: Y=1#¥ 
P=. TRUE. 
X(N)}EN#ZETA (ALPHA#M/PI)}+(1.,-1.) 


Table 2 indicates which type of expression may be equated to each type of variable in an arithmetic state- 
ment. D indicates that the assignment is performed directly (no conversion of any sort is done); R indicates 
that only the real part of the variable is set to the value of the expression (the imaginary part is set to zero);. 
C means that the expression is converted to the type of the variable; and H means that only the high-order 


portion of evaluated expression is assigned to the variable. 


The expression value is made to agree in type with the assignment variable before replacement occurs. 


For example, in the statement: 
THETA=W«(ABETA+E) 


if THETA is an integer and the expression is real, the expression value is truncated to an integer before 


assignment to THETA. 


TABLE 2 ALLOWED ASSIGNMENT STATEMENTS 


Expression 


Logical, 
: D t 
Variable Complex ouble Oc! al, or 
Precision Literal 
Constant 


i. ae D,H,L C,H,L R,D,H,L D,H,L 
recision 


D - Direct Replacement 
C - Conversion between integer and floating point 


R - Real only 


| - Set imaginary part to 0 
H - High order only 


L - Set low order part to 0 


CHAPTER 4 


CONTROL STATEMENTS 


FORTRAN compiled programs normally execute statements sequentially in the order in which they were 
presented to the compiler. However, the following control statements are available to alter the normal 
sequence of statement execution: GOTO, IF, DO, PAUSE, STOP, END, CALL, RETURN. CALL and 


RETURN are used to enter and return from subroutines. 


GO TO STATEMENT 


The GO TO statement has three forms: unconditional, computed, and assigned. 


Unconditional GO TO Statements 
Unconditional GO TO statements are of the form: 
GOTOn 


where n is the number of an executable statement. Control is transferred to the statement numbered n. 


Computed GO TO Statements 
Computed GO TO statements have the form: 
GO TO (Ny Mo, a “My dei 
where NypMoreee eM are statement numbers, and i is an integer expression. 
This statement transfers control to the statement numbered My eMor ee ee My if it has the value 1, 2, ...,k, 


respectively. If i exceeds the size of the list, or is zero, execution will proceed to the next executable 


statement. 


For example, in the statement: 
GO TO (20,10, 5),K 


the variable K acts as a switch, causing a transfer to statement 20 if K=1, to statement 10 if K=2, or to 


statement 5 if K=3. 


Assigned GO TO Statement 
Assigned GO TO statements have two forms: 


GOTOk 
and 


GO TO k(n) MorM3, Sie) 


where k is a nonsubscripted integer variable and NyrMor++ +n) are statement numbers. Both forms of the 
assigned GO TO have the effect of transferring control to the statement whose number is currently asso- 
ciated with the variable k. This association is established through the use of the ASSIGN statement, the 


general form of which is: 


ASSIGN i TO k 


If more than one ASSIGN statement refers to the same integer variable name, the value assigned by the 


last executed statement is the current value. 
EXAMPLES: ASSIGN 21 TO INT ASSIGN 1000 TO INT 


GO TO INT GO TO INT, (2,21, 1000, 310) 


IF STATEMENT 


IF statements have two forms in FORTRAN IV: numerical and logical . 


Numerical IF Statements 
Numerical IF statements are of the form: 
IF (expression) yoy 


where nye n are statement numbers. 


o°"3 
This statement transfers control to the statement numbered Aye Nor Ng if the value of the numeric expression 


is less than, equal to, or greater than zero, respectively. The expression may not be complex. 


EXAMPLES: IF (ETA)4,7,12 
IF (KAPPA-L(10))20, 14,14 


Logical IF Statements 


Logical IF statements have the form: 
IF (expressions)S 


where S is a complete statement. 


The expression must be logical. S may be any executable statement other than a DO statement or another 


logical IF statement. 
If the value of the expression is .FALSE., control passes to the next sequential statement. 


If the value of the expression is . TRUE., statement S is executed. After execution of S, control passes 
to the next sequential statement unless S is a numerical IF statement or a GO TO statement; in these 


cases, control is transferred as indicated. 


If the expression is . TRUE. and S is a CALL statement, control is transferred to the next sequential state- 


ment upon return from the subroutine. 


Numbers are present in the logical expression: 


IF (B)Y=X«SIN(Z) 
W=Y ««2 


If the value of B is .TRUE., the statements Y=X«*SIN(Z) and W=Y««2 are executed in that order. If the 
value of B is .FALSE., the statement Y=X*SIN(Z) is not executed. 


EXAMPLES: IF (T.OR.S)X=Y+1 
IF (Z.GT.X(K)) CALL SWITCH (S,Y) 
IF (K.EQ.INDEX) GO TO 15 


NOTE: Care should be taken in testing floating point numbers for equality in 
IF statements as rounding errors may cause unexpected results. 


DO STATEMENT 
The DO statement simplifies the coding of iterative procedures. DO statements are of the form: 


DO n i=m),m 


a3 


where n is a statement number, i is a nonsubscripted integer variable, and m),mo,m., are any integer 


expressions. If m, is not specified, it is understood to be |. 


3 


a. ee 


The DO statement causes the statements which follow, up to and including the statement numbered n, to 
be executed repeatedly. This group of statements is called the range of the DO statement. The integer 
variable i of the DO statement is called the index. The values of My, Mo, and.m., are called, respectively, 


the initial, limit, and increment values of the index. 


A zero increment (m,)is notallowed. The increment may be negative if m >m,- Ifm)<mo, the increment 


m,, must be positive. The parameters m, and m, may have like or unlike signs as long as ms is always 


3 


larger than m3, Or mM, is always larger than mo: 


EXAMPLES: Form Restriction 
DO 10 IF1,5,2 
DO 10 I=5,1,-1 
DO 10 IFJ,K,5 J<K 
DO 10 l=J,K,-5 J>K 
DO 10 L=1,J,-K I< J,K<0 or I>J,K>0 
DO 10 L=I,J,K I<J,K>0 or I>J,K<O 


Initially, the statements of the range are executed with the initial value assigned to the index. This 
initial execution is always performed, regardless of the values of the limit and increment. After each 
execution of the range, the increment value is added to the value of the index and the result is compared 
with the limit value. If the value of the index is not greater than the limit value, the range is executed 
again using the new value of the index. When the increment value is negative, another execution will 


be performed if the new value of the index is not less than the limit value. . 


After the last execution of the range, control passes to the statement immediately following the range. 


This exit from the range is called the normal exit. Exit may also be accomplished by a transfer from within 
the range. The value of the index after normal exit is indeterminate 
since the compiler stores and increments the index in index register 
15 whenever possible. DDT users should note that they may have to 
examine register 15 at a breakpoint inside a DO loop to determine the 
current value of the index. 

The range of a DO statement may include other DO statements, provided that the range of each contained 


DO statement is entirely within the range of the containing DO statement. That is, the ranges of two DO 


statements must intersect completely or not at all. A transfer into the range of a DO statement from out- 


side the range is not allowed. 


Within the range of a DO statement, the index is available for use as an ordinary variable. After a trans- 
fer from within the range, the index retains its current value and is available for use as a variable. The 


values of the initial, limit, and increment variables for the index and the value of the index 
itself, may not be altered within the range of the DO statement. 


The range of a DO statement must not end with a GO TO type statement or a numerical IF statement. A 
logical IF statement is allowed as the last statement of the range. In this case, control is transferred as 


follows. The range is considered ended when, and if, control would normally pass to the statement follow- 
ing the entire logical IF statement. 


As an example, consider the sequences: 


DO 5 K=1,4 
5 IF(X(K). GT. ¥(K))Y¥(K}=X(K) 
en : 


Statement 5 is executed four times whether the statement Y(K}X(K) is executed or not. 
Statement 6 is not executed until statement 5 has been executed four times. 


EXAMPLES: DO 22 L=1,30 
DO 45 K=2,LIMIT,-3 
DO 7 X=T,MAX,L 


CONTINUE STATEMENT 


The CONTINUE statement has the form: 

CONTINUE 
This statement is a dummy statement, used primarily as a target for transfers, particularly as the last state- 
ment in the range of a DO statement. For example, in the sequence: 


DO 7 K=START, END 
IF (X(K))22, 13,7 
7 CONTINUE 


a positive value of X(K) begins another execution of the range. The CONTINUE provides a target address 


for the IF statement and ends the range of the DO statement. 


PAUSE STATEMENT 


The PAUSE statement enables the program to incorporate operator activity into the sequence of automatic 


events. The PAUSE statement assumes.one of three forms: 


PAUSE 
PAUSE n 
PAUSE 'xxxxx' 


where n is an unsigned string of six or less octal digits, and 'xxxxx' is a literal message. 


Execution of the PAUSE statement causes the message or the octal digits, if any, to be typed on the 


user's teletypewriter. Program execution may be resumed (at the next executable FORTRAN statement) 


from the console by typing "G," followed by a carriage return. Program execution may be terminated by 


typing "X," followed by carriage return. 


EXAMPLE: PAUSE 167 
PAUSE 'NOW IS THE TIME' 


STOP STATEMENT 
The STOP statement has the form: 


STOP 


The STOP statement terminates the program and returns control to the monitor system. (Termination of a 


program may also be accomplished by a CALL to the EXIT or DUMP subroutines. ) 


END STATEMENT 
The END statement has the form: 


END 


The END statement informs the compiler to terminate compilation and must be the physically last statement 


of the program. 


CHAPTER 5 


INPUT/OUTPUT STATEMENTS AND SUBROUTINES ‘eo 


Input/output statements are used to control the transfer of data between computer memory and peripheral 


devices and to specify the format of the output data. Input/output statements may be divided into three 


categories, as follows: 


1. Nonexecutoble statements that enable conversions between internal form data 
within core memory and external form data (FORMAT), or specify lists of arrays and | 


variables for input/output transfer (NAMELIST). These statements are compiled ° 


"in line" with jumps around the body of the statement. Therefore,, 


attempts to execute these statements will simply result in transfer 
to the next statement. 


2. Statements that specify transmission of data between computer memory and I/O 


devices: READ, WRITE, PRINT, PUNCH, TYPE, ACCEPT. 


3. Statements that control magnetic tape unit mechanisms: REWIND, BACKSPACE, 
END FILE, UNLOAD, SKIP RECORD. 


4. Statements and subroutines that control directory devices: 


REWIND, BACKSPACE, END FILE, SKIP RECORD, CALL IFILE, CALL OFILE, 
CALL RREAD, CALL RWRI. 


NONEXECUTABLE STATEMENTS 


The FORMAT statement enables the user to specify the form and arrangement of data on the selected ex- 
ternal medium. The NAMELIST statement provides for conversion and input/output transmission of data 


without reference to a FORMAT statement. 


FORMAT Statement 


FORMAT statements may be used with any appropriate input/output medium. FORMAT statements are of » 
the form: 


1 


1 1 * 
n FORMAT (Spp Sor ++ +S /Sq 1 Soy S/o. .) 


where n is a statement number, and each S is a data field specification. 


FORMAT statements may be placed anywhere in the source program. Unless the FORMAT statement con- 
tains only alphanumeric data for direct input, output transmission, it will be used in conjunction with the . 


list of a data transmission statement. . Q 


Slashes are used to specify unit records . 


There are no restrictions on record lengths. The programmer 


is responsible for restricting output line length for devices 


which have a fixed line length. 


During transmission of data, the object program scans the designated FORMAT statement. If a specification 
for a numeric field is present (see "Input/Output Lists" p 35) and the data transmission statement contains 
items remaining to be transmitted; transmission takes place according to the specification. This process 
ceases and execution of the data transmission statement is terminated as soon as all specified items have 
been transmitted. Thus, the FORMAT statement may contain specifications for more items than are speci- 
fied by the data transmission statement. Conversely, the FORMAT statement may contain specifications 


for fewer items than are specified by the data transmission statement. 


The following types of field specifications may appear in a FORMAT statement: numeric, numeric with 
scale factors, logical, alphanumeric. The FORMAT statement also provides for handling multiple record 
formats, formats stored as data, carriage control, skipping characters, blank insertion, and repetition. 
If an input list requires more characters than the input device supplies for a given unit record, blanks 


are supplied. 


Numeric Fields 


Numeric field specification codes and the corresponding internal and external forms of the numbers are 


listed in Table 3. 
The conversions of Table 3 are specified by the forms: 


1}. Dw.d 

Ew.d 

Fw.d 

Iw 

. Ow 

Gw.d (for real) 


Oo WW & W KN 


Gw (for integer or logical) 


Gw.d,Gw.d (for complex) 


respectively. The letter D, E, F, 1, O, or G designates the conversion type; w is an integer specifying 


the field width, which may be greater than required to provide for blank columns between numbers; d is 


we OO oS 


an integer specifying the number of decimal places to the right of the decimal point or, for G conversion, 
the number of significant digits. (For D, E, F, and G input, the position of the decimal point in the 


external field takes precedence over the value of d in the format.) 


For example, 
FORMAT (15, F10.2,D18.10) 
could be used to output the line, 
bbb32bbbb—17.60bbb . 5962547681 D+03 


on the output listing. 


The field width w should always be large enough to include spaces for the decimal point, sign, and expo- 
nent. In all numeric field conversions if w is not large enough to accommodate the converted number, the 
excess digits on the left will be lost; if the number is less than w spaces in length, the number is right- 


adjusted in the field. 


TABLE 3 NUMERIC FIELD CODES 


Conversion 


Code External Form 


D Binary floating point Decimal floating point 
double-precision with D exponent 

E Binary floating point Decimal floating point 

with E exponent 

F Binary floating point Decimal fixed point 

| Binary integer Decimal integer 
Binary integer Octal Integer 
One of the following: Single precision 
single precision decimal floating point, 
binary floating point, integer, logical (T or 
binary integer, F), or complex (two 
binary logical, or decimal floating point 
binary complex numbers), depending 


upon the internal form 


Numeric Fields with Scale Factors 
Scale factors may be specified for D, E, F, and G conversions. A scale factor is written nP where P is 


the identifying character and n is a signed or unsigned integer that specifies the scale factor. 


For F type conversions (or G type, if the external field is decimal fixed point), the scale factor specifies 


a power of ten so that 


external number = (internal number)* jo Keele factor) 


For D, E, and G (external field not decimal fixed point) conversions, the scale factor multiplies the 
number by a power of ten, but the exponent is changed accordingly leaving the number unchanged except 


in form. For example, if the statement: 
FORMAT (F8.3,E16.5) 
corresponds to the line 
bb26 .451 bbbb—0.41321E—01 
then the statement 
FORMAT (—1PF8.3, 2PE16.5) 
might correspond to the line | 


bbb2.645bbb—41 .32157E-03 


In input operations, F type (and G type, if the external field is decimal fixed point) conversions are the 


only types affected by scale factors. 


When no scale factor is specified, it is understood to be zero. However, once a scale factor is specified, 
it holds for all subsequent D, E, F, and G type conversions within the same format unless another scale 
factor is encountered. The scale factor is reset to zero by specifying a scale factor of zero. Scale factors 


have no effect on | and O type conversions. 


Logical Fields 
Logical data can be transmitted in a manner similar to numeric data by use of the specification: 


Lw 


where L is the control character and w is am integ” specifying the field width. The data is transmitted 


as the value of a logical variable in the input/output list. 


If on input, the first nonblank character in the data field is T or F, the value of the logical variable will 


be stored as true or false, respectively. If the entire data field is blank, a value of false will be stored. 


On output, w minus | blanks followed by T or F will be output if the value of the logical variable is true 


or false, respectively. 


Variable Field Width 


The D, E, F, G, I, and O conversion types may appear without the specification of the field width w. 


In the case of input, omitting the w implies that the numeric field is delimited by any character which 


would otherwise be illegal in the field in addition to the characters blank and carriage return 


Provided they follow the numeric field. For example, input according to 


the f ct: 
orma 10 FORMAT(21,F,E, ©) 


might appear as: 
-10,3,15.621,.0016E-10,777 


On output, omitting the w has the following effect: 


Format Becomes 
D D25.16 
E E15.7 
F F15.7 
G G15.7 or G25.16 
I 115 
O O15 


Alphanumeric Fields 

Alphanumeric data can be transmitted in a manner similar to numeric data by use of the form Aw, where 
A is the control character and w is the number of characters in the field. The alphanumeric characters 
are transmitted as the value of a variable in an input/output list. The variable may be of any type. For 
the sequence: 


READ 5,V 
5 FORMAT (A4) 


causes four characters to be read and placed in memory as the value of the variable V. 


Although w may have any value, the number of characters transmitted is limited by the maximum number 
of characters which can be stored in the space allotted for the variable. This maximum depends upon the 
variable type. For an integer or real variable the maximum is five characters; for a double precision or 
complex variable, the maximum is ten characters. If w exceeds the maximum, the leftmost characters are 
lost on input and replaced with blanks on output. If, on input, w is less than the maximum, blanks are 
filled in to the right of the given characters until the maximum is reached. If, on output, w is less than 


the maximum, the leftmost w characters are transmitted to the external medium. 


Alphanumeric Data Within Format Statements 
Alphanumeric data may be transmitted directly into or from the format statement by two different methods: 


H-conversion, or the use of single quotes. 


In H-conversion, the alphanumeric string is specified by the form nH. H is the control character and n is 
the number of characters in the string counting blanks. For example, the format in the statement below 


can be used to print PROGRAM COMPLETE on the output listing. 


FORMAT (17H PROGRAM COMPLETE) 


Referring to this format in a READ statement would cause the 17 characters to be replaced with a new 


string of characters. 


The same effect is achieved by merely enclosing the alphanumeric data in quotes. The result is the same 
as in H-conversion; on input, the characters between the quotes are replaced by input characters, and, 
on output, the characters between the quotes (including blanks) are written as part of the output data. A 


quote character within the data is represented by two successive quote marks. For example, referring to: 
FORMAT ('DON''T') 


with an output statement would cause DON'T to be printed. 


Mixed Fields 

An alphanumeric format field may be placed among other fields of the format. For example, the statement: 
FORMAT (15, 7H FORCE=F10.5) 

can be used to output the line: 


bbb22bFORCE=bb1 7.68901 


The separating comma may be omitted after an alphanumeric format field, as shown above. 


Complex Fields 


Complex quantities are transmitted as two independent real quantities. The format specification consists 


of two successive real specifications or one repeated real specification. For instance, the statement: 
FORMAT (2E15.4, 2(F8.3,F8.5)) 


could be used in the transmission of three complex quantities. 


Repetition of Field Specifications 


Repetition of a field specification may be specified by preceding the control character D, E, F, 1, O, G, 


L, or A by an unsigned integer giving the number of repetitions desired. For example: 
FORMAT (2E12.4, 315) 
is equivalent to: 


FORMAT (E12.4,E12.4, 15,15, 15) 


Repetition of Groups 


A group of field specifications may be repeated by enclosing the group in parentheses and preceding the 


whole with the repetition number. For example: 
FORMAT (218, 2(E15.5, 2F8.3)) 
is equivalent to: 


FORMAT (218, £15.5, 2F8.3, £15.5, 2F8.3) 


Multiple Record Formats 


To handle a group of input/output records where different records have different field specifications, a 


slash is used to indicate a new record. For example, the statement: 
FORMAT (308/15, 2F8.4) 
is equivalent to: 
FORMAT (308) 
for the first record and 


FORMAT (15, 2F8. 4) 


for the second record. 


The separating comma may be omitted when a slash is used. When n slashes appear at the end or beginning 
of a format, n blank records may be written on output or records skipped on input. When n slashes appear 


in the middle of a format, n-1 blank records are written or n-1 records skipped. 


Both the slash and the closing parenthesis at the end of the format indicate the termination of a record. 
If the list of an input/output statement dictates that transmission of data is to continue after the closing 


parenthesis of the format is reached, the tormat is repeated from the last open parenthesis of level one or 


zero. Thus, the statement: 


FORMAT (F7.2, (2(E15.5,€15.4),17)) 


level 0 level J level J level 0 


causes the format: 
F7.2,2(E15.5,E15.4),17 
to be used on the first record, and the format: 
2(E15.5,E15.4), 7 


to be used on succeeding records. 


As a further example, consider the statement: 


FORMAT (F7.2/(2(E15.5,E15.4), 17) ) 


The first record has the format: 
F7.2 


and successive records have the format: 


2(E15.5,E15.4),17 


Formats Stored as Data 


The ASCII character string comprising a format specification may be stored as the values of an array. 
Input/output statements may refer to the format by giving the array name, rather than the statement number 
of a FORMAT statement. The stored format has the same form as a FORMAT statement excluding the word 


"EORMAT." The enclosing parentheses are included. 


As an example, consider the sequence: 


DIMENSION ARRAY (2) 

READ 1, (ARRAY(I), I=1,2) 
FORMAT (2A4) 

READ ARRAY,K,X 


The first READ statement enters an ASCII string into the array ARRAY. . In the second READ statement, 


ARRAY js referred to as the format governing conversion of K and X. » 
Carriage Control 5 


The first character of each ASCII record controls the spacing of the line printer or Teletype. This character 
is usually set by beginning a FORMAT statement for an ASCII record with IHa, where a is the desired con- 


trol character. The line spacing actions, listed below, occur before printing: 


Character Effect 
space skip to next line 
0 skip a line 2 
form feed - go to top of next page, This is simulated 


On user consoles which do not have hardware 
form feeds. 


+ suppress skipping - will overprint line 
- skip 2 lines 


Characters not in the table will act as spaces. However, characters 
will be added to the table from time to time. 


A $ (dollar sign) as a format field specification code Suppresses the 
carriage-return at the end of the line. 


Spacing 


Input and output can be made to skip forward to any position within a 
FORTRAN record by use of the format code: 


Tw 


where T is the control character and w is an unsigned integer constant specifying the position in a FORTRAN 


record where the transfer of data is to begin. 


For example, 


2FORMAT(T3@,'BLACK'T59@,'WHITE') 


would cause the following line to be printed: 


Print Position 3@ Print Position 5@ 
J J 
BLACK WHITE 


For input, the statements: 
1 FORMAT(T35, 'MONTH') 
READ (3,1) 


would cause the first 34 characters of the input data to be skipped, and the next 5 characters would re- 


place the characters M, O, N, T, and H in storage. 


Blank or Skip Fields 


Blanks may be introduced into an output record or characters skipped on an input record by use of the 
specification nX. The control character is X; n is the number of blanks or characters skipped and must 


be greater than zero. For example, the statement: 
FORMAT (5H STEPI5, 1OX2HY=F7.3) 


may be used to output the line: 


bSTEPbbb28bbbbbbbbbbY=b—3. 872 


NAMELIST Statement 


The NAMELIST statement, when used in conjunction with special forms of the READ and WRITE statements, 
provides a method for transmitting and converting data without using a FORMAT statement or an I/O list. 


The NAMELIST statement has the form: 
NAMELIST/X,/A, Ao, me -A/X/B, Bo, sas -,B.. . SK LC, Co, ne Co 


where the X's are NAMELIST names, and the A's, B's, and C's are variable or array names. 


Each list or variable mentioned in the NAMELIST statement is given the NAMELIST name immediately 
preceding the list. Thereafter, an I/O statement may refer to an entire list by mentioning its NAMELIST 


name. For example: 


NAMELIST/FRED/A, B, C/MARTHA/D,E 


states that A, B, and C belong to the NAMELIST name FRED, and D and E belong to MARTHA. 
The use of NAMELIST statements must obey the following rules: 


1. A NAMELIST name may not be longer than six characters; it must start with an 
alphabetic character; it must be enclosed in slashes; it must precede the list of entries 


to which it refers; and it must be unique within the program. 


2. A NAMELIST name may be defined only once and must be defined by a NAMELIST 
statement. After a NAMELIST name has been defined, it may only appear in READ or 
WRITE statements. The NAMELIST name must be defined in advance of the READ or 
WRITE statement. 


3. A variable used in a NAMELIST statement cannot be used as a dummy argument in 


a subroutine definition. 


4. Any dimensioned variable contained in NAMELIST statement must have been de- 


fined in a DIMENSION statement preceding the NAMELIST statement. 


Input Data for NAMELIST Statements 


When a READ statement refers to a NAMELIST name, the first character of all input records is ignored. 
Records are searched until one is found with a $ or & as the second character immediately followed by the 
NAMELIST name specified. Data is then converted and placed in memory until the end of a data group is 
signaled by a $ or & either in the same record as the NAMELIST name, or in any succeeding record as 
long as the $ or & is the second character of the record. Data items must be separated by commas and be 


of the following form: 


V=Ky,Ko,-- eK 


where V may be a variable name or an array name, with or without subscripts. The K's are constants 
which may be integer, real, double precision, complex (written as (A, B) where A and B are real), or 
logical (written as T or .TRUE., and F or .FALSE). A series of J identical constants may be represented 
by J*K where J is an unsigned integer and K is the repeated constant 'ogical and complex constants 
must be equated to logical and complex variables, respectively. The other types of constants (real, double 
precision, and integers) may be equated to any other type of variable (except logical or complex), and 
will be converted to the variable type. For example, assume A is a two-dimensional real array, B is a 


one-dimensional integer array, C is an integer variable, and that the input data is as follows: 


$FRED A(7, 2 4,B=3,6%2.8, C=3.32$ 

Column 2 
A READ statement referring to the NAMELIST name FRED will result in the following: the integer 4 will 
be converted to floating point and placed in A(7,2). The integer 3 will be converted to floating point 


and placed in B(1) and the floating point number 2.8 will be placed in B(2), B(3),..., B(7). The floating 


point number 3.32 w be converted to the integer 3 and placed in C. 


Output Data for NAMELIST Statements — 


When a WRITE statement refers to a NAMELIST name, all variables and arrays and their values belonging 
to the NAMELIST name will be written out, each according to its type. The complete orray is written 
out by columns. The output data will be written so that: 


1. The fields for the data will be large enough to contain all the significant digits. 


2. The output can be read by an input statement referencing the NAMELIST name. 


For example, if JOE is a 2x3 array, the statements: 


NAMELIST/NAM1/JOE, K1, ALPHA 
WRITE (u, NAM1) 


will generate the following form of output: 


Column 2 
$NAME! 
JOE= -6.75, . 234E-04, 68.0, 
-17.8, 0.0, -~.197E+07, 
K1=73.1, ALPHA=3 


DATA TRANSMISSION STATEMENTS 


The data transmission statements accomplish input/output transfer of data that may be listed in a NAMELIST 
statement or defined in a FORMAT statement. When a FORMAT statement is used to specify formats, the 
data transmission statement must contain a list of the quantities to be transmitted. The data appears on the 


external media in the form of records. 


Input/Output L ists 


The list of an input/output statement specifies the order of transmission of the variable values. During 
input, the new values of listed variables may be used in subscript or control expressions for variables 


appearir iater in the list. For example: 
READ 13,L, ACL) ,B(L+1) 
reads a new value of L and uses this value in the subscripts of A and 38. 


The transmission of array variables may be controlled by indexing similar to that used in the DO statement. 


The list of controlled variables, followed by the index control, is enclosed in parentheses. For example, 


READ 7, (X(K),K=1,4),A 


is equivalent to: 


READ 7,X(1),X(2),X(3),X(4),A 


As in the DO statement, the initial, limit, and increment values may be given as integer expressions: 


READ 5, N, (GAIN(K), K=1,M/2,N) 


The indexing may be compounded as in the following: e 


READ 11, ( (MASS(K,L), K=1,4),L=1, 5) 


The above statement reads in the elements of array MASS in the following order: 
MASS(1,1), MASS(2,1),..., MASS(4,1), MASS(1,2),...,MASS(4, 5) 
If an entire array is to be transmitted, the indexing may be omitted and only the array identifier written. 


The array is transmitted in order of increasing subscripts with the first subscript varying most rapidly. Thus, 


the example above could have been written: 


READ 11, MASS 


Entire arrays may also be designated for transmission by referring to a NAMELIST name (see description of e 


NAMELIST statement). 


Input/Output Records 


All information appearing on external media is grouped into records. The maximum amount of information 


in one record and the manner of separation between records depends upon the medium. For punched cards, 


each card constitutes one record; on a teletypewriter a record is one line, and so forth. The amount of 
information contained in each ASCII record is specified by the FORMAT reference and the I/O list. For 
& magnetic tape binary records, the amount of information is specified by the I/O list. 


Each execution of an input or output statement initiates the transmission of a new data record. Thus, the 


statement: 
READ 2, FIRST,SECOND, THIRD 
is not necessarily equivalent to the statements: 


READ 2, FIRST 
, READ 2, SECOND 
READ 2, THIRD 


* since, in the second case, at least three separate records are required, whereas, the single statement: 
READ 2, FIRST, SECOND, THIRD 
may require one, two, three, or more records depending upon FORMAT 2. 
If an input/output statement requests less than a full record of information, the unrequested part of the 


record is lost and cannot be recovered by another input/output statement without repositioning the record. 


& If an input/output list requires more than one ASCII record of information, successive records are read. 


PRINT Statement 


PRINT refers to a logical device named "PRINT." If a physical device 
is not assigned the logical name "PRINT" before execution, the error 
message "DEVICE PRINT NOT AVAILABLE" will be given by the FORTRAN 
operating system. There are two forms for the print statement: 

PRINT f, list 


PRINT f 


: where f is a format reference. 


The data is converted from internal to external form according to the designated format. If the data to be 


transmitted is contained in the specified FORMAT statement, the second form of the statement is used. 


EXAMPLES: PRINT 16,T, (B(K), K=1,M) 
PRINT F106, SPEED, MISS 


In the second example, the format is stored in array F106. 


PUNCH Statement 


PUNCH refers to a logical device named "PUNCH. If a physical device 
is not assigned the logical name "PUNCH" before execution, the error 
message "DEVICE PUNCH NOT AVAILABLE" will be given by the FORTRAN 


Operating System. There are two forms for the PUNCH statements: 
PUNCH f, list 
PUNCH f 


where f is a format reference. 


Conversion from internal to external data forms is specified by the format reference. If the data to be 


transmitted is contained in the designated FORMAT statement, the second form of the statement is used. 


EXAMPLES: PUNCH 12,A,B(A),C(B(A) ) 
PUNCH 7 


TYPE Statement 
The TYPE statement assumes one of two forms: 


TYPE f, list 
TYPE f 


where f is a format reference. 


This statement causes the values of the variables in the list to be read from memory and listed on the user's 
teletypewriter. The data is converted from internal to external form according to the designated format. 


If the data to be transmitted is contained in the designated FORMAT statement, the second form of the 


statement is used: 


EXAMPLES: TYPE 14,K, (A(L), L=1, K) 
TYPE FMT 


WRITE Statement 
The WRITE statement assumes one of the following forms: 
WRITE(u, f) list 
WRITE (u, f) 
WRITE(u, N) 
WRITE (u) list 


where u is a unit designation, f is a format reference, and N is a NAMELIST name. 


The first form of the WRITE statement causes the values of the variables in the list to be read from memory 
and written on the unit designated in ASCII form. The data is converted to external form as specified by 


the designated FORMAT statement. 


The second form of the WRITE statement causes information to be read directly from the specified format 


and written on the unit designated in ASCII form. 


The third form of the WRITE statement causes the names and values of all variables and arrays belonging 
to the NAMELIST name, N, to be read from memory and written on the unit designated. The data is 


converted to external form according to the type of each variable and array. 


The fourth form of the WRITE statement causes the values of the variables in the list to be read from mem- 


ory and written on the unit designated in binary form. 


READ Statement 


The READ statement assumes one of the following forms: ~ 


READ f, lis . These forms require 

READ f that a physical device be 
assigned the logical name 

READ\u, f) list "CARD" before execution. 

READ (u, f) 

READ(u, N) 

READ(u) list 


where f is a format reference, u is a unit designation, and N is a NAMELIST name. 


The first form of the READ statement causes information to be read from cards and put in memory as values 
of the variables in the list. The data is converted from external to internal form as specified by the ref- 


erenced FORMAT statement. 


EXAMPLE: READ 28, Z1,Z2,Z3 


The second form of the READ statement is used if the data read from cards is to be transmitted directly into 


the specified format. 


EXAMPLE: READ 10 


The third form of the READ statement causes ASCII information to be read from the unit designated and 


stored in memory as values of the variables in the list. The data is converted to internal form as specified 


by the referenced FORMAT statement. 
EXAMPLE: READ(1,15)ETA, PI 


The fourth form of the READ statement causes ASCII information to be read from the unit designated and 


transmitted directly into the specified format. 
EXAMPLE: READ(N, 105) 


The fifth form of the READ statement causes data of the form described in the discussion of input data for 
NAMELIST statements to be read from the unit designated and stored in memory as values of the variables 


or arrays specified. 
EXAMPLE: READ(2, FRED) 


The sixth form of the READ statement causes binary information to be read from the unit designated and 


stored in memory as values of the variables in the list. 


EXAMPLE: READ(M)GAIN, Z, Al 


ACCEPT Statement 
The ACCEPT statement assumes one of two forms: 


ACCEPT f, list 
ACCEPT f 


where f is a format reference. 
This statement causes information to be input from the user's teletypewriter and put in memory as values 


of the variables in the list. The data is converted to internal form as specified by the format. If the 


transmission of data is directly into the designated format, the second form of the statement is used. 


EXAMPLES: ACCEPT 12, ALPHA, BETA 
ACCEPT 27 


DEVICE CONTROL STATEMENTS AND SUBROUTINES 


TABLE 4 DEVICE CONTROL STATEMENTS 


Statement Effect 


BACKSPACE unit . Backspaces one logical record for either 
ASCII or binary files. A single back- 
space is legal for any device that does 
input. Repeated backspaces are only 
permitted on directory devices and mag- 
netic tapes. Backspace is a read 
operation only. 


END FILE unit Closes the current file for input and 
output. If output is active, a physical 
end-of-file record is written for mag- 
netic tape. For directory devices, 

END FILE is identical to REWIND. 


REWIND unit © Does an END FILE for any device and then 
rewinds if the device is a magnetic tape. 


SKIP RECORD unit Skips one logical record for ASCII or 
binary files. 


UNLOAD unit Same as REWIND except that if the device 
is a magnetic tape, it will unload. 


TABLE 4A DEVICE CONTROL SUBROUTINES 


Statement Effect 


IFILE (unit, name, ext) Opens file "NAME.EXT" for input on a 


directory device. Name is an ASCII 
variable consisting of one to five alpha- 
numeric characters. EXT is an ASCII 


variable consisting of one to three char- 
acters. EXT is an optional argument with 
DAT as the default argument. 


OFILE (unit, name, ext) Opens a new file "NAME.EXT" for output 
on a directory device. Arguments are 
the same as for IFILE. 


RREAD (unit, record, count) Positions the read pointer so that the > | 
next READ statement will read record 
number "RECORD" each of which are of 
length “COUNT”. This is only appli- 
cable for binary files of constant 
record length. (See Appendix 7) 


RWRI (unit, record, count) Converse of RREAD, i.e. positions write 
pointer so that the next WRITE state- 
ment will write record number "RECORD" 
of length "COUNT". The file can only 
be extended one record at a time, e.g. P 
if the file is of length 20 records, 
record 21 can be written, but record 2? 
cannot. (See Appendix 7) 


The following calls set or clear flags defining conditions for the 
file associated with unit "UNIT". These calls are made available 
to allow the programmer to have maximum control over the I/O. 


SETBIN (unit) Sets a flag defining file to be in 
binary format. This flag is automati- 
cally set by a binary read or write 
operation. It is only needed if an 
operation such as SKIP RECORD is the 
first operation on a binary file. 


SETASC (unit) Sets a flag defining the file to be 
ASCII. Since any file is assumed to 
be ASCII until a binary operation is 
done, the only application of this 
call would be in switching modes 
within a file. 


SETRAN (unit) Sets a flag defining the file to be a 
random access binary file of constant 
record length. This flag is automati- 
cally set by a random read or write 
operation. It is only needed if an 
operation such as SKIP RECORD is the “ 
first operation ona file. 


CLRRAN (unit) This clears the random file flag. . 
SETEFT (unit) Sets flag for end-of-file testing. 

This flag is also set by a call to 

EOFC. 


CLREFT (unit) Clears end-of-file test flag. 
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SETSEQ (unit) 


CLRSEQ (unit) 


TABLE 4B 


Statement 


EOFC (unit) 


Sets flag to allow sequence numbers to 
be read as ASCII data. A sequence num- 
ber consists of six characters (five 
digits followed by a tab). (See Appen- 
dix 7) 


Clears flaz and causes sequence numbers 
to be skipped on input. 


DEVICE CONTROL FUNCTIONS 


Effect 


Returns a true (-1) [false (@)] if an 
end-of-file has [has not] been read on 
unit “UNIT". EOFC calls SETEFT. 

(See Appendix 7) 


CHAPTER 6 


SPECIFICATION STATEMENTS 


Specification statements allocate storage and furnish information about variables and constants to the 


compiler. Specification statements may be divided into three categories, as follows: 


1. Storage specification statements: DIMENSION, COMMON, and EQUIVALENCE. 
2. Data specification statements: DATA and BLOCK DATA. 


3. Type declaration statements: INTEGER, REAL, DOUBLE PRECISION, COMPLEX, 
LOGICAL, SUBSCRIPT INTEGER, and IMPLICIT. 


The following specification statements, if used, appear in the program prior to any executable statement: 


DIMENSION statement 

EXTERNAL statementt 

COMMON statement 

EQUIVALENCE statement 

Type declaration statements 

Arithmetic function definition statements! 
DATA statements 

IMPLICIT statements 


In addition, arrays must be dimensioned before being referenced in a NAMELIST or DATA statement. 
STORAGE SPECIFICATION STATEMENTS 


DIMENSION Statement 


The DIMENSION statement is used to declare identifiers to be array identifiers and to specify the number 
and bounds of the array subscripts. The information supplied in a DIMENSION statement is required for 


the allocation of memory for arrays. Any number of arrays may be declared ina single DIMENSION 


statement. The DIMENSION statement has the form: 


DIMENSION Sy , S57 Stag 5. 


where S is an array specification. 


> EXTERNAL and arithmetic function definition statements are described in Chapter 7. 


Each array variable appearing in the program must represent an element of an array declared in a 
DIMENSION statement, unless the dimension information is given ina COMMON or TYPE statement. 


Dimension information may appear only once for a given variable. 


Each array specification gives the array identifier and the minimum and maximum values which each of its 


subscripts may assume in the following form: 


identifier(min/max, min/max, ...,min/max) 
The minima and maxima must be integers. The minimum must not exceed the maximum. 


For example, the statement: 
DIMENSION EDGE(-1/1, 4/8) 


specifies EDGE to be a two-dimensional array whose first subscript may vary from —] to 1 inclusive, and 


the second from 4 to 8 inclusive. 


Minimum values of 1 may be omitted. For example, 
NET(5, 10) 
is interpreted as: 


NET(1/5, 1/10) 


EXAMPLES: DIMENSION FORCE(-1/1,0/3,2,2,-7/3) | 
DIMENSION PLACE(3, 3, 3), JI(2, 2/4), K(256) 


Arrays may also be declared in the COMMON or type declaration statements in the same way: 


COMMON X(10,4),Y,Z 
INTEGER A(7, 32),B 
DOUBLE PRECISION K(-2/6, 10) 


Adjustable Dimensions 


Within either a FUNCTION or SUBROUTINE subprogram, DIMENSION statements may use integer vari- 
ables in an array specification, provided that the array name and variable dimensions are dummy arguments 
of the subprogram. The actual array name and values for the dummy variables are given by the calling 
program when the subprogram is called. The variable dimensions may not be altered within the subprogram 


and must be less than or equal to the explicit dimensions declared in the calling program. 


EXAMPLE: SUBROUTINE SBR(ARRAY, M1 ,M2,M3, M4) 
DIMENSION ARRAY (M1/M2,M3/M4) 


DO 27 L=M3,M4 
DO 27 K=M1,M2 


27 ARRAY (K, LF VALUE 


END 


The calling program for SBR might be: 


DIMENSION AI(10, 20), A2(1000, 4) 
CALL SBR(AI,5, 10, 10, 20) 


CALL SBR(A2, 100, 250, 2, 4) 


END 


ALLOCATING ARRAYS AT EXECUTION TIME 


Arrays that are not stored in common can be defined at execution 

time. This enables a programmer to vary the size of his core image: 

as the size of his problem changes, instead of having to allow for 

the largest possible problem. In addition, space used by the loader 

at lead time can be recovered and used if the user is near the limit 

for core. To use this feature, the user writes his main program as 

a subroutine and uses variable dimensioning for the arrays to be de- . 
fined at execution. A new main program is written which inputs the 

array sizes from the console and calls the subroutine "ALLOT". The 

call to “ALLOT" is: 


CALL ALLOT(NAME,N1,N2,....,-1,P1,P2,....P(I-1)) 
where name is the subroutine name of the old main program, the "N's" 


are array sizes, "I" is the number of additional arguments (includ- 
ing itself), and the "P's" are additional parameters to be passed. 


Example: 


SUBROUTINE OLDM(ARAY1,ARAY2,ARAY3,K1,K2,K3) 
DIMENSION ARAY1(K1),ARAY2(K2,K2) ,ARAY3(K2,K3) 


oe? 


The new main program to call ALLOT might be as follows: 


EXTERNAL OLDM 

ACCEPT 900,K1,K2,K3 

N1=K1 

N2=K2*K2 

N3=K2*K3 

CALL ALLOT(OLDM,N1,N2,N3,-4,K1,K2,K3) 
900  FORMAT(6L) 

END 


When the new main program is run, ALLOT will calculate the start- 
ing address of each array from free storage and pass control to 
the program "OLDM". 


COMMON Statement 


The COMMON statement causes specified variables or arrays to be stored in an area available to other 
programs. By means of COMMON statements, the data of a main program and/or the data of its subpro- 


grams may share a common storage area. 


The common area may be divided into separate blocks which are identified by block names. A block is 


specified as follows: 

/block identifier/identifier, identifier, ..., identifier 
The identifier enclosed in slashes is the block name. The identifiers which follow are the names of the 
variables or arrays assigned to the bloc and are placed in the block in the order in which they appear 


in the block specification. A common block may have the same name as a variable in the same program 


or as any subroutine or function name in the same job. 


The COMMON statement has the general form: 


COMMONS BLOC KI/A, B, C/BLOCK2/D,E,F/... 


where BLOCK], BLOCK2, ... are the block names, and A,B,C, ... are the variables to be 


assigned to each block. For example, the statement: 
COMMON/R/X, Y,1T/C/U,V,W,Z 


indicates that the elements X,Y, and T are to be placed in block R in that order, and that U,V,W, and Z 


are to be placed in block C. 
Block entries are linked sequentially throughout the program, beginning with the firsts COMMON statement. 
For example, the statements: 


COMMON/D/ALPHA/R/A, B/C/S 
COMMON/C/X, Y/R/U, V,W 


have the same effect as the statement: 


COMMON/D/ALPHA/R/A, B,U,V,W/C/S,X,Y 


One block of common storage, referred to as blank common, may be left unlabeled. Blank common is 


indicated by two consecutive slashes. For example, 


COMMON/R/X, Y//B,C,D 


indicates that B, C, and D are placed in blank common. The slashes may be omitted when blank common 


is the first block of the statement: 


COMMON B,C,D 


Storage allocation for blocks of the same name begins at the same location for all programs executed to- 


gether. For example, if a program contains 
COMMON A,B/R/X,Y,Z 
as its firsts COMMON statement, and a subprogram has 
COMMON/®/U,V,W//D,E, F 


as its first COMMON statement, the quantities represented by X and U are stored in the same location. 


A similar correspondence holds for A and D in blank common. 


Common blocks may be any length provided that no program attempts to enlarge a given common block 


declared by a previously loaded program. 


Array names appearing in COMMON statements may have dimension information appended if the arrays 


are not declared in DIMENSION or type declaration statements. For example, 
COMMON ALPHA, T(15, 10,5), GAMMA 
specifies the dimensions of the array T while entering T in blank common. Variable dimension array 


identifiers may not appear in a COMMON statement, nor may other dummy identifiers. 


Each array name appearing in a COMMON statement must be dimensioned somewhere in the program con- 


taining the COMMON statement. 


EQUIVALENCE Statement 


The EQUIVALENCE statement causes more than one variable within a given program to share the same 


storage location. The EQUIVALENCE statement has the form: 
EQUIVALENCE(V, Vor nF (Vere Pe Pee 


where the V's are variable names. 


The inclusion of two or more references in a parenthetical list indicates that the quantities in the list are 


to share the same memory location. For example, 
EQUIVALENCE (RED, BLUE) 


specifies that the variables RED and BLUE are stored in the same place. 


The relation of equivalence is transitive; e.g., the two statements, 


EQUIVALENCE(A, B), (B,C) 
EQUIVALENCE(A,8,C) 


have the same effect. 
The subscripts of array variables must be integer constants. 


EXAMPLE: EQUIVALENCE (X, A(3), Y(2, 1, 4)), (BETA(2, 2), ALPHA) 


EQUIVALENCE and COMMON 


identifiers may appear in both COMMON and EQUIVALENCE statements provided the following rules are 


observed. 


1. No two quantities in common may be set equivalent to one another. 


uh a a 


2. Quantities placed in a common block by means of EQUIVALENCE statements may 


cause the end of the common block to be extended. 


For example, the statements: 


COMMON/R/X, Y,Z 
DIMENSION A(4) 
EQUIVALENCE(A, Y) 


causes the common block R to extend from X to A(4), arranged as follows: 


xX 
Y A(l) (same location) 
Z A(2) (same location) 
A(3) 
A(4) 


3. EQUIVALENCE statements which cause extension of the start of a common block 


are not allowed. For example, the sequence: 


COMMON/R/X, Y, Z 
DIMENSION A(4) 
EQUIVALENCE(X, A(3) ) 


is not permitted, since it would require A(1) and A(2) to extend the starting location 


of block R. 


DATA SPECIFICATION STATEMENTS 


The DATA statement is used to specify initial or constant values for variables. The specified values are 
compiled into the object program, and become the values assumed by the variables when program execution 


begins. 


DATA Statement 


The data to be compiled into the object program is specified in a DATA statement. The DATA statement 


has the form: 
DATA list/d, ido, bas J Nist/d dy ye Maree: 


where each list is in the same form as an input/output list, and the d's are data items for each list. 


Indexing may be used ina list provided the initial, limit, and increment (if any) are given as constants. 


Expressions used as subscripts must have the form: 


c,xitc 


1 2 


where Cy and cy are integer constants and i is the induction variable. If an entire array is to be defined, 
only the array identifier need be listed. Variables in common may appea or the lists only if the DATA 


statement occurs in a BLOCK DATA subprogram. 


The data items following each list correspond one-to-one with the variables of the list. Each item of the 


data specifies the value given to its corresponding variable. 


Data items may be numeric constants, alphanumeric strings, octal constants, or logical constants. For 


example, 


DATA ALPHA, BETA/5, 16.£-2/ 


specifies the value 5 for ALPHA and the value .16 for BETA. 


Alphanumeric data is packed into words according to the data word size in the manner of A conversion; 
however, excess characters are not permitted. The specification is written as nH followed by n characters 


or is imbedded in single quotes. 


Octal data is specified by the letter O or the character ", followed by a signed or unsigned octal integer 


of one to twelve digits. 
Logical constants are written as .TRUE.,.FALSE., T, or F. 
EXAMPLE: DATA NOTE, K/7HRADIANS, O-7712/ 


Any item of the data may be preceded by an integer followed by an asterisk. The integer indicates the 


number of times the item is to be repeated. For example: 
DATA(A(K), K=1, 20)/61E2, 19*32E1/ 


specifies 20 values for the array A; the value 6100 for A(1); the value 320 for A(2) through A(20). 


BLOCK DATA SUBPROCRAM 


The BLOCK DATA statement has the form: 


BLOCK DATA 


This statement declares the program which follows to be a data specification subprogram. Data may be 


entered into common only. 


The first ssatement of the subprogram must be the BLOCK DATA statement. The subprogram may contain 


only the declarative statements associated with the duta being defined. 


EXAMPLE: BLOCK DATA 
COMMON/?/S, Y/C/Z,W,V 
DIMENSION Y(3) 

COMPLEX Z 


DATA Y/1E-1, 2*3E2/, X, Z/11.877D0, (—1.41421,1.41421)/ 
END 


Data may be entered into more than one block of common in one subprogram. 


TYPE DECLARATION STATEMENTS 


The type declaration statements INTEGER, REAL, DOUBLE PRECISION, COMPLEX, LOGICAL, IMPLICIT, 
and SUBSCRIPT INTEGER are used to specify the type of identifiers appearing in a program. An identifier 


may appear in only one type statement. Type statements may be used to give dimension specifications 
for arrays. 


The explicit type declaration statements have the general form: 


type identifier, identifier, identifier... 
where type is one of the following: 
INTEGER, REAL, DOUBLE PRECISION, COMPLEX,LOGICAL, 
SUBSCRIPT INTEGER 


The listed identifiers are declared by the statement to be of the stated type. Fixed-point variables ina 
SUBSCRIPT INTEGER statement must fall between -27~” and 2°”. 


IMPLICIT Statement 


The IMPLICIT statement has the form: 


IMPLICIT type, (a, rns SN ars types (az,a4, | 


where type represents one of the following: INTEGER, REAL, LOGICAL, COMPLEX, DOUBLE PRECISION; 


and a,0,, ... represent single alphabetic characters, each separated by commas, or a range of characters 


(in alphabetic sequence) denoted by the first and last characters of the range separated by a minus sign 
(e.g., (A-D)). 


This statement causes any program variable which is not mentioned in a type statement, and whose first 
character is one of those listed, to be typed according to the type appearing before the list in which the 


character appears. As an example, the statement: 
IMPLICIT REAL(A—D,L, N—P) 


causes all variables starting with the letters A through D, L, and N through P to be typed as real, unless 


they are explicitly declared otherwise. 


The initial state of the compiler is set as if the statement 
IMPLICIT REAL(A—H, O-—Z), INTEGER(I-N) 
were at the beginning of the program. This state is in effect unless an IMPLICIT statement changes the 
above interpretation; i.e., identifiers, whose types are not explicitly declared, are typed as follows: 
1. Identifiers beginning with I, J, K, L, M, or N are assigned integer type. 
2. Identifiers not assigned integer type are assigned real type. 
If the program contains an IMPLICIT statement, this statement will override throughout the program the 


implicit state initially set by the compiler. No program may contain more than one IMPLICIT declaration 


for the same letter. 
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CHAPTER 7 


SUBPROGRAM STATEMENTS 


FORTRAN subprograms may be either internal or external. Internal subprograms are defined and may be 
used only within the program containing the definition. The arithmetic function definition statement is 


used to define internal functions. 


External subprograms are defined separately from (i.e., external to) the programs that call them, and are 
complete programs which conform to all the rules of FORTRAN programs. They are compiled as closed 
subroutines; i.e., they appear only once in the object program regardless of the number of times they arv 


used. External subprograms are defined by means of the statements FUNCTION and SUBROUTINE. 


Dummy Identifiers 


Subprogram definition statements contain dummy identifiers, representing the arguments of the subp:ogram. 
They are used as ordinary identifiers within the subprogram definition and indicate the sort of arguments 
that may appear and how the arguments are used. The dummy identifiers are replaced by the actual argu- 


ments when the subprogram is executed. 


Library Subprograms 


The standard FORTRAN IV library for the PDP-10 includes built-in functions, FUNCTION subprograms, 
and SUBROUTINE subprograms, listed and described in Appendixes 1, 2, and 3, respectively. Built-in 
functions are open subroutines; that is, they are incorporated into the object program each time they are 
referred to by the source program. FUNCTION and SUBROUTINE subprograms are closed subroutines; 


their names derive from the types of subprogram statements used to define them. 


ARITHMETIC FUNCTION DEFINITION STATEMENT 
The arithmetic function definition statement has the form: 


identifier (identifier, identifier, ...)- expression 


This statement defines an internal subprogram. The entire definition is contained in the single statement. 


The first identifier is the name of the subprogram being defined. 


Arithmetic function subprograms are single-valued functions with at least one argument. The type of the 


function is determined by the type of the function identifier. 


The identifiers enclosed in parentheses represent the arguments of the function. These are dummy identifiers; 
they may appear only as scalar variables in the defining expression. Dummy identifiers have meaning and 
must be unique only within the defining statement. Dummy identifiers must agree in order, number, and 


type with the actual arguments given at execution time. 


Identifiers, appearing in the defining expression, which do not represent arguments are treated as ordinary 
variables. The defining expression may include external functions or other previously defined arithmetic 


statement functions. . 
All arithmetic function definition statements must precede the first executable statement of the program. 


EXAMPLES: SSQR(K}K«(K+1 )#(2eK+1)/6 
ACOSH (X}(EXP(X/A}t+ EXP(—X/A) )/2 


In the last example above, X is a dummy identifier and A is an ordinary identifier. At execution time, the 


function is evaluated using the current value of the quantity represented by A. 
FUNCTION SUBPROGRAMS 


A FUNCTION subprogram is a single-valued function that may be called by using its name as a function 
name in an arithmetic expression, such as FUNC(N), where FUNC is the name of the subprogram that eval- 
vates the corresponding function of the argument N. A FUNCTION subprogram begins with a FUNCTION 
statement and ends with an end statement. It returns control to the calling program by means of one or 


more RETURN statements. 


FUNCTION Statement 


The FUNCTION statement has the form: 


FUNCTION identifer(argument, argument, .. .) 


This statement declares the program which follows to be a function subprogram. The identifier is the name 
of the function being defined. This identifier must appear as a scalar variable and be assigned a value dur- 


ing execution of the subprogram which is the function value. 


Arguments appearing in the list enclosed in parentheses are dummy arguments representing the function 
argument. The arguments must agree in number, order, and type with the actual arguments used in the 


calling program. Function subprogram arguments may be expressions, alphanumeric strings, array names, 


or subprogram names. 


Dummy arguments may appear in the subprogram as scalar identifiers, array identifiers, or subprogram 
identifiers. A function must have at least one dummy argument. Dummy arguments representing array names 
must appear within the Gbercgran ina DIMENSION statement, or one of the type statements that provide 
dimension information. Dimensions given as constants must equal the dimensions of the corresponding arrays 
in the calling program. In a DIMENSION statement, dummy identifiers may be used to specify adjustable 


dimensions for array name arguments. For example, in the statement sequence: 


FUNCTION TABLE(A,M, N,B,X, Y) 
_ DIMENSION A(M,N), B(10), C(50) 


The dimensions of array A are specified by the dummies M and N, while the dimension of array B is given 
as a constant. The various values given for M and N by the calling program must be those of the actual 
arrays which the dummy A represents. The arrays may each be of different size but must have two dimen- 


sions. The arrays are dimensioned in the programs that use the function. 


Dummy dimensions may be given only for dummy arrays. In the example above the array C must be given 
absolute dimensions, since C is not a dummy identifier. A dummy identifer may not appear in an EQUIV- 


ALENCE statement in the function subprogram . 


- A function must not modify any arguments which appear in the FORTRAN arithmetic expression calling the 
function. Modification of implicit arguments from the calling program, such as variables in common and 
DO loop indexes, is not allowed. The only FORTRAN statements not allowed in a function subprogram are 


SUBROUTINE, BLOCK DATA, and another FUNCTION statement. 


Function Type 


The type of the function is the type of identifier used to name the function. This identifier may be typed, 
implicitly or explicitly, in the same way as any other identifier. Alternatively, the function may be ex- 
plicitly typed in the FUNCTION statement itself by replacing the word FUNCTION with one of the 


following: 


INTEGER FUNCTION 

REAL FUNCTION 

COMPLEX FUNCTION 
LOGICAL FUNCTION 

DOUBLE PRECISION FUNCTION 


Y For example, the statement: 
COMPLEX FUNCTION HPRIME(S, N) 
is equivalent to the statements: 


FUNCTION HPRIME(S, N) 
COMPLEX HPRIME 


EXAMPLES: | FUNCTION MAY(RANGE, EP, YP, ZP) 
| COMPLEX FUNCTION COT(ARG) 
- DOUBLE PRECISION FUNCTION LIMIT(X, Y) 


SUBROUTINE SUBPROGRAMS 


A subroutine subprogram may be multivalued and can be referred to only by a CALL statement. A sub- 
routine subprogram begins with a SUBROUTINE statement and returns control to the calling program by 


means of one or more RETURN statements. 


SUBROUTINE Statement 
& The SUBROUTINE statement has the form: 


SUBROUTINE identifier(argument, argument, .. .) 


This statement declares the program which follows to be a subroutine subprogram. The first identifier is 
the subroutine name. The arguments in the list enclosed in parentheses are dummy arguments representing 
the arguments of the subprogram. The dummy arguments must agree in number, order, and type with the 


actual arguments used by the calling program. 


Subroutine subprograms may have expressions, alphanumeric strings, array names, and subprogram names 


as arguments. The dummy arguments may appear as scalar, array, or subprogram identifiers. 


Dummy identifiers which represent array names must be dimensioned within the subprogram by a DIMENSION 
or type declaration statement. As in the case of a function subprogram, either constants or dummy identi- 
fiers may be used to specify dimensions in a DIMENSION statement. The dummy arguments must not appear 


in an EQUIVALENCE or COMMON statement in the subroutine subprogram. 


A subroutine subprogram may use one or more of its dummy identifiers to represent results. The subprogram 


name is not used for the return of results. A subroutine subprogram need not have any argument at all. 


EXAMPLES: SUBROUTINE FACTOR(COEFF,N,ROOTS) 
SUBROUTINE RESIDU(NUM, N,DEN,M,RES) 
SUBROUTINE SERIES 


The only FORTRAN statements not allowed in a function subprogram are FUNCTION, BLOCK DATA, 
and another SUBROUTINE statement. 


CALL Statement 
The CALL statement assumes one of two forms: 


CALL identifier 


CALL identifier (argument ,argument, ... ,argument) 


The CALL statement is used to transfer control to subroutine subprogram. The identifier is the subprogram 


name. 


The arguments may be expressions, array identifiers, alphanumeric strings or subprogram identifiers; argu- 
ments may be of any type, but must agree in number, order, type, and array size (except for adjustable 
arrays, as discussed under the DIMENSION statement) with the corresponding arguments in the SUBROUTINE 


statement of the called subroutine. Unlike a function, a subroutine may produce more than one value and 


cannot be referred to as a basic element in an expression. 


A subroutine may use one or more of its arguments to return results to the calling program. If no arguments 


at all are required, the first form is used. 


EXAMPLES: CALL EXIT 
CALL SWITCH(SIN, 2. LE.BETA, X««4, Y) 
CALL TEST(VALUE, 123, 275) 


The identifier used to name the subroutine is not assigned a type and has no relation to the types of the 


arguments. Arguments which are constants or formed as expressions must not be modified by the subroutine. 


RETURN Statement 
The RETURN statement has the form: 
RETURN 


This statement returns control from a subprogram to the calling program. Normally, the last statement 


executed in a subprogram is a RETURN statement. Any number of RETURN statements may appear in a © | 


subprogram. 


EXTERNAL Statement 


Function and subrou' ‘ne subprogram names may be used as the actual arguments of subprograms. Such 
subprogram na. .es must be distinguished from ordinary variables by their appearance in an EXTERNAL state- 


ment. The EX~CRNAL statement has the form: 


EXTERNAL identifier, identifier, ..., identifier 


This statement declares the listed identifiers to be subprogram names. Any subprogram name given as an 


argument to another subprogram must appear in an external declaration in the calling program. 


EXAMPLE: EXTERNAL SIN,COS 
CALL TRIGF(SIN, 1.5, ANSWER) 
CALL TRIGF(COS, .87, ANSWER) 


END 


SUBROUTINE TRIGF(FUNC, ARG, ANSWER) 
ANSWER = FUNC(ARG) 


RETURN 
END 


To reference external variables from a MACRO-10 program, place the variables in named COMMON. 
Use the name of the variable as the name of the COMMON block: - 


COMMON /A/A, /B/B (13),/C/, C(6,7) 


General Form Page References 
ASSIGN 4OM  <cbueeambeleaedyoia swash ert waede ka anneeey meee 17 
CALL name (0) ,G5,-26) cee eeee eee ese eee eeeeceteneeeereeeceenes 56 
CONTINGE  siuiswundyt edi cote ialnanre en eueewenen unica 20 
DO i MAMMA pM eee eee e eee eeeeeenn eee eeeeeeeeeneneeeeeas 18 
OCOWOL, ~psca densa eae cave asuae cone sewoeeoewe ess ehae seca 16 
COMO!  ssdgrageisnstan miaadnmucestenerseesuueeeenyaseds 17 
COTOMALGAIOE udorascancuamcsriecintitenaneeneons 17 
GO TO (i -i5, pM: adeeeeeecs Siw nese iene Gh Re data a ate bw Ova Butte 16 
IFC )iprignig cece eee eee eee ee eee e eee e eee e eee e eee eeeens 17 
HES «aise tece natn ees eo einstein 18 
PAUSE aah ben ead sh Cammew eee cee Ue haere sau wena ume iads 20 
PAUSE ciiieentapeuiiegeutcneet we taauaaeee amet eu eenes 21 
PROSE | ins lisivasutaltien hte co cua oltia w eahigcnloaaieelaiw aly eet arena esata 21 
RETURN: Ha ecete ak antes eiSedecsa te Seeanuhnreaaueee amie 21 
STOR Cs ae uGewhtwe eared aca peewee te Negras somocmeees 21 
END! gadivit anu ogee teat Baesas Ueraaaea es eae eA eS 21 

INPUT/OUTPUT STATEMENTS 
General Form | Page References 
ACCERIC | dtcuotatesetneertudadeedueed Mambueweeaaameae eee. 38 
PCCERD EMS ars tesa nee evans aiatid eatea aba alechansceemlen 38 
BACKSPACE URI | sist cle bagi tite Sow deren cota ease 39 
END: FILE Umit) owt eiete rete reek ove edaw et ekaneaaitad ews 39 
PORMAIG)” _-sc4 Deel pee oua men wah ee Gu ieee eee a 22 
PRUNES: (suv, gue wane nmne cin seak ee ante eG ae xed ea ne ewe baie eae ae 35 
PRUNE AMIS okt Sake en ea eee ann ede o ae andes baane sree Se hae ae 35 
FUNGHE:  siws ueuieuudececwamo vetted we oleue nel sia akese abe 36 


APPENDIX 1 


SUMMARY OF ALC FORTRAN IV STATEMENTS 


CONTROL STATEMENTS 


General Form Page Reference 


READE . coccuacenauaceieuaieinee acetate amo ateak Mneoanades 37 
READ f, list ict Riise a unceueaae wa enna Vesamaeans es teeter aaa 37 
READ (unit, f) see we ee eeeeceeece icnaiab at shave tneaesedeas 37 
READ (URIF,ONIS  Xeamuntesdectcancuseaelodeahotearenctaenees a 37 
READ (unit)list ws eee eee irate lack Sees Wes waco ss Pi fecaucecg 2 37 
READ (unit, name ,) Seen Te ee re eT re eee ee eer eee foes 37 
REWIND unit" Ceo a eee tone ete iee Te 
SKIPRECORD Unit sc ce dinciveessdeusessaucsareaees soewbetan 39 
FVPED. .cntuciieee cc decec banter aeietsdoadnaeasuiwea ves ion 36 
TYPES isk, <a usiweidssedcaonateaves seh ie keds e ame senna 36 
WRITE QUGIEY  teitncas ude tctautis gard Gauloenaiaatvieatuuonaenes 36 
WRITEAUDIL ENE: ireundoteciti wae cca sehu seatuneseduaswnonad 36 
WRITERIDIISE Gas Sa hieern avis woh ie se edness ashen 36 
WRITE (unit, name ;) bce oe aaes bande ete a LO eau ioe Dew Ase ene ee 36 
UNLOAD unit eee ee eee eee eee 39 


INPUT/OUTPUT SUBROUTINES 


General Form . Page Reference 


IFILE (unit, name, ext) 
OFILE (unit, name, ext) 


eoeoeseeeeereeeeoe eee eee eee eanee eee ee.e 39 


eeeoseeteeeteeeereeseeoeeeeeeenwneeenvneeeene 39 


RREAD (unit, record, count) oenteoeeeeeweeseee eee eeeee eee 40 
RWRI (unit, record, count) eoeeeteeeeeeeeeeseeeeeeeeee eee 40 
SETBIN (unit), oeeoeoec0seseeetveeseeeeeeeeereeeeeeeepeeeesee eoee35u8eee@¢°@ 40 
SETASC (unit), oeeeo3e#eeenreereeeeeeeeeeevsxsesvsee eoeeoeeeevs#e54#:e#@¢ @ 40 
SETRAN (unit), oeeoeoeeeseesvse @ eee eeene eeeeeetkeeeee@ oes @ eee te 40 
CLRRAN (unit) eeee#e#e#e#8e e@eee7eevseecs3540vle*eeseesveeeeeseeeseeeeee @eeteee 40 
SETEFT (unit), eoeeese0e8seteeee3ee3erseeeefeeeeteeeweeeeeee oe eoeeeseve 40 
CLREFT (unit), oee0ees3e«#steeteenrteeeeeereeeeeeeteesteeeeee eoee35sc5eeee¢ee 40 
SETSEQ (unit) ooo eee e ewe eee ewe ememceml mC OCC OCH HO PO HO Hh HOO Oo eho ov ee#e?e 41 
CLRSEQ (unit) eoeeese5seeer#rke8eee oeoeoeeoewere eee eee wm woe ee wm wm wee 41 
INPUT/OUTPUT FUNCTIONS 
General Form Page Reference 
EOFC (unit) eoeeveeseeeeeeeees “oe eee eee ee ee we eee eo wore e®e¢nseeene 41 


SPECIFICATION STATEMENTS 


General Form 
_ TE I ea i a aE 


COMMON a(ny Mo, -+-)eb(ng,n4, eee), oe 
COMPLEX A(ny Mor ++ )eb(ngsn yy ee s)y oe 
DATA tUr ee Ski rkorkg, «+ / 

VeWr ee Sk ark kyy eee/e. 
DIMENSION Sa A 


DOUBLE PRECISION a(n, Mores -)bing,n4, ee pee 


EQUIVALENCE a(n, pe -), bro, eda) fein, ppesivn 
(c(n,, -+),d(ny, Arey rey Ree 

EXTERNAL y,z,... 

IMPLICIT type, (1 \7 15), types (14-14), iets 

INTEGER a(n) ,n, ee -)b(ng ny, CS ee 

LOGICAL a(n, More -)bing,n4, odd voted 

NAMELIST /name,/a,b, ts -/name,/c,d, Sieg 

REAL a(n, Mores )birgsn4, ee re 

SUBSCRIPT INTEGER a(n, Mores -),b(n,, Aes) Dene 
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Page Reference 
45 
50 
48 


42 
30 
47 


57 
50 
50 
50 
31 


50 
50 
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ARITHMETIC STATEMENT FUNCTION DEFINITION 


General Form Page Reference 
name(a,b,...)=e 51 
NOTE: 
a, 1Oo, es are expressions 
a,b,c,d are variable names 
e . is an expression 
e) is a noncomplex expression 
e, is a logical expression 
f is a format number 
g , is a format specification 
"h' is an alphanumeric string 
i,i, ‘ ios es are statement numbers 
j is an integer constant 
ky rko, awe are constants of the general form jxk where 
k is any constant 
l los ue are letters | 
list is an input/output list 
m is an integer variable name 
mi ,M5,M. are integer expressions © 
Ny Moroes are dimension specifications 
name is a subroutine or function name 
name, ,name., are NAMELIST names 
s is a statement (not DO or logical IF) 
t,U,V,W are variable names or input/output lists 
type,, types, diets are type specifications 
unit is an integer variable or constant specifying a 
logical device number 
yz - are external subprogram names 


APPENDIX 2 


FORTRAN IV LIBRARY FUNCTIONS 


This appendix contains descriptions of all standard function subprograms provided with the FORTRAN IV 
library for the PDP-10. These functions may be called by using the function mnemonic as a function name 


in an arithmetic expression. 


Function Definition Number of Misi Type of | 
Arguments Argument Function 


Absolute value | Arg| ] ABS Real Real 
Absolute value | Arg| 1 [ABS Integer Integer 
Absolute value | Arg | ! DABS Double Double 
Truncation Sign of Arg times ] AINT Real Real 
largest integer <|Arg| 
Truncation Sign of Arg times 1 INT Real Integer 
largest integer <|Arg| 
Truncation Sign of Arg times l IDINT Double Integer 
largest integer <|Arg| 
Remaindering Arg (mod Argo) t 2 AMOD Real Real 
Remaindering Arg (mod Arg.) t 2 MOD Integer Integer 
Choosing largest Max(Arg, Ardy: sey >2 AMAXO Integer Real 
value ; AMAX 1 Real Real 
MAXO Integer Integer 
MAX 1 Real Integer 
DMAX | Double Double 
Choosing Min(Arg, Ao, eso) 22 AMINO Integer Real 
smallest value AMIN] Real Real 
MINO Integer Integer 
MIN} Real Integer 


DMIN1 Double Double 


The function MOD or AMOD (a1, a9) is defined as a; - la)/a9] a5, where [x] is the integer whose 
magnitude does not exceed the magnitude of x and whose sign is the same as x. 


5 | ; Number of Type of 
Transfer of sign Sgn(Arg.)* | sonldralara| =| SIGNN SIGN Real Real 
Transfer of sign Sgn(Argo)s | Arg, | ISIGN Integer Integer 
Transfer of sign Sgn(Argo)s | Arg, | DSIGN Double Double 
Positive Arg, -Min(Arg, Argo) DIM Real Real 
difference | 
Positive Arg, -Min(Arg, sArgo) IDIM Integer Integer 
difference 
Complex For Arg=XtiY , C=X-iY CONJG Complex Complex 
conjugate 
Conversion from FLOAT Integer Real 
integer to real 
Conversion from Result is largest FIX Real Integer 
real to integer integer <a 
Express single . DBLE Real Double 
precision argu- 
ment in double 
precision form, 
low order part =0 
Express two real C=Arg, tixArg, CMPLX Real Complex 
arguments in 
complex form 
Obtain most SNGL Double Real 
significant part . 
of double pre- 
cision argument 
Obtain real REAL Complex Real 
part of complex 
argument 
Obtain AIMAG Complex Real 
imaginary part 
of complex 
argument 

; Arg 
Exponential e EXP Real Real 


Function 


Natural 
logarithm 


Common 
logarithm 


Arc-sine 
Arc-cosine 
Arctangent 
Arctangent of 
the quotient 
of two 
arguments 
Sine (radians) 


Sine (degrees) 


Cosine 
(radians) 


Cosine 
(degrees) 


Hyperbolic 
tangent 


Hyperbolic 


sine 


Hyperbolic 
cosine 


Square root 


Remaindering t 


Exponential 


Natural 
logarithm 


Common 
logarithm 


log, o(Ara) 


asin(Arg) 
acos(Arg) 
atan(Arg) 


atan(Arg,/Arg,) 


sin(Arg) 
sin(Arg) 


cos(Arg) 
cos(Arg) 
tanh(Arg) 
sinh(Arg) 
cosh(Arg) 


(Arg) 7 2 


Arg, (mod Argo) 
ag 


log. (Arg) 


log) (Arg) 


1 ALOG 

1 ~  ALOGIO- 

i ASIN 

1 ACOS 

1 ATAN 

2 ATAN2 

1 SIN 

1 SIND 

1 cos 

1 COSD 

1 TANH 

1 SINH 

1 COSH 

1 —SQRT- 

2 DMOD 

1 DEXP 

1 DLOG 
DLOGIO 


Real 
Real 


Real 
Real 
Real 


Real 


Real 
Real 


Real 
Real 
Real 
Real 
Real 


Real 


Double 


Double 


Double 


Double 


Argument 


Type of 
Function 


Real 
Real 


Real 
Real 
Real 
Real 
Real 


Real 


Double 


Double 


Double 


Double 


+The function DMOD (a1, 99) is defined as a, - [a)/a9] a2, where [x] is the integer whose magnitude 
does not exceed the magnitude of x and whose sign is the same as the sign of x. 


i nis 


Function 


Arctangent 


Arctangent 
of two 
arguments 


Sine(radians) 


Cosine 
(radians) 


Square root 
Absolute value 


Exponential 


Natural 
logarithm 


Complex sine 


Complex 
cosine 


Complex 
square root 


Arguments 


atan(Arg) 


atan(Arg, /Argy) 


sin(Arg) 


cos(Arg) 


(Arg) '/2 


c=ox2ey4)/2 


Arg 


e 
log, (Arg) 
sin(Arg) 


cos(Arg) 


caxsiy)/? 
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DATAN2 


DSIN 


DCOS 


DSQRT 


CABS 


CEXP 


CLOG 


CSIN 


-CCOS 


CSQRT 


Type of 
Argument — Function 
Double Double 
Double Double 
Double Double 
Double Double 
Double Double 
Complex Real 
Complex Complex 
Complex Complex 
Complex Complex 
Complex Complex 
Complex Complex 


APPENDIX 3 


FORTRAN IV LIBRARY SUBROUTINES 


This appendix contains descriptions of all standard subroutine subprograms provided within the FORTRAN IV 


library for the PDP-10. These subprograms are closed subroutines and may be called with a CALL statement. 


Subroutine Name Effect 
EXIT Returns control to the monitor and, therefore, terminates the . 
execution of the program. 
% 
Subroutine Name Effect 
SLITE(i) Where i is an integer expression, turns sense lights on or off. 
For 1<i<36 sense light i will be turned on. If i=0, all sense 
lights will be turned off. 
SLITET(i, j) Checks the status of sense light i and sets the variable j ac- | 


cordingly and turns off sense light i. If iis on, j is set to 1; 


and if i is off, j-is set to 2. 


OVERFL(j) Checks the status of the AR OV flag and sets the variable j 
accordingly. If the AR OV flag is on, j is set to 1. If the 


flag is off, j is set to 2. 


APPENDIX 4 


ALC FORTRAN IV OPERATING SYSTEM 


SUBPROGRAM CALLING SEQUENCES 


FORTRAN Subroutines 
FORTRAN subroutine calling sequences appear as follows: 


JSA 16, NAME 
ARG CODE, 1A, 


ARG CODE, ,A, 


where NAME is the name of the subroutine, ARG is a "pseudo-op" equivalent to a JUMP instruction 
(a "no-op"), CODE, CODE,,, etc. are 4-bit codes with the values: 


0 integer argument 
] Unused 
2 Real argument 
3 Logical argument 
4 Octal argument 
5 Literal argument 
6 Double precision argument 
7 Complex argument 
and A, Ao. etc. are the argument addresses. 


All accumulators are saved in subroutines except 0 for subroutines with single-word arguments and except 
0 and | for subroutines with double-word arguments (high order or real part in 0 and low order or imaginary 
part in 1). All scalar arguments in a subroutine call are transferred into and restored from the subroutine 


by value. 


FORTRAN Function Subprograms and Library Functions 
The FORTRAN function calling sequence is the same as that for subroutines. The function value is returned 
in accumulator 0 or accumulators 0 and 1. Scalar arguments are not restored as their values may not be 


modified within a function. 


INPUT/OUTPUT 


The ALC FORTRAN IV library contains a fortran operating system 
(FORSE.) which controls fortran input/output. Communication be- 
tween the program and the operating system is accomplished sy a 
number of defined operation codes which are listed in Tale 5. 
These codes can be used by a macro programmer or by a :ortran 
programmer wishing to make patches to his program. 


TABLE 5 DEFINED OPERATION CODES 


Operation Code Name Function 
(OCTAL) 
15 RESET. Reinitializes all 1/0. First instruc- 
tion in a fortran program. 
16 IN. Set up to read ASCII records. 
17 OUT. Set up to write ASCII records. 
20 DATA. Transfer ASCII information in or out. 
21 FIN. Finish binary or ASCII read or write. 
22 RTB. Read binary records. 
23 WTB. Write binary records. 
24 MTOP. Device control statements. 
25 SLIST. Short list I/0. (ARRAY 1/0). 
26 INF. Open file for input. 
27 OUTF. Open file for output. 
31 NLI. Name list input. 
32 NLO. Name list output. 
33 RRB. Read random access binary file. 
34 WRB. Write random access binary file. 
35 GIO. Get I/O channel and buffer header. 


FORTRAN UNL7T - I/0 DEVICE ASSIGNMENTS 


Fortran unit numbers coricspond to a logical or physical device 
according to a table in the operating system named DEVTB. (Fig.1) 
Device TTY refers to the user's console. FILE "N" logical de- 
vices correspond to file drawers in the ALC directory systen. 
FILE "N" devices are special in that a file name must be asso- 
ciated with each one before it can be used. To explicitly assign 
file names to file drawers the IFILE or OFILE statement is used. 
If an explicit assignmert is not made, a default file name of 
FILE "N".DAT is used. 


The other logical devices are named for the physical devices they 
represent: READ (card reader), PUNUn (card punch), PRINT (line 
printer). These devices can be used by assigning the physical de- 
vice the proper logical name. 


The monitor command: 

-ASSIGN LPsr PRINT 

would cause a subsequent fortran PRINT statement to do output on 
the line printer. Assignments can also be used to move output 
around. If fortran unit three were used for output and it were 
desirable to put the output temporarily on the user's console, 
the monitor command 


.ASSIGN TTY FILE3 


would do it. 


FIGURE 1 DEVICE TABLE 


Unit Device Comment 
1 FILE1 Directory system file drawers. 
2 FILE2 ; 
3 FILE3 
4 FILE4 
5 TTY User console 
6 TTY User console 
7 FILE7 
18 FILE18 
FORTRAN STATEMENT DEVICE 
PRINT PRINT 
PUNCH PUNCH 
TYPE TTY 
ACCEPT TTY 
READ READ 
ASCII mode 


On output, blocks are packed for every device. The block size 
varies from 16 for the user console to 128 for directory devices 
and magnetic tapes. The first character of every logical record 

is a format control character which is interpreted for line 
printers and user consoles and written out directly for other 
devices. To list ASCII files which have been written with a for- 
tran program, the format characters must be interpreted at the time 
of listing. Line printers and user consoles also cause output at 
the end of every WRITE statement. 
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On input, the carriage return character or altmode (ESC) is the 
record mark. 


RECORD FORMATING 
BINARY MODE 


In binary mode records are packed in 128 word blocks on magnetic 
tape and on directory devices. Record length is arbitrary. The 
first word of each record (or piece of a record if the record 

contains more than 126 data words) is a control word of the form: 


_ W N 


Where "W'' is the word count of the data words in this block and 
"N" is @ in all but the last block, in which case "N" is the num- 
ber of blocks in the record. (A record contains all the data 
corresponding to one READ or WRITE statement.) 


A constant overhead of one word per record (the control word) 
plus one word per block (extra control word for split records) 
is assured by inserting zero words when the two overhead words 
coincide. This is to facilitate the implementation of random 
access to binary files with fixed record length. Therefore, an 
array of size 126 would fit exactly into one physical block. 


REREAD 


For every READ statement there is a corresponding REREAD statement 
which will read the last record again. 


BUFFERING 


All I/0 is double buffered including random access I/O. The sys- 
tem is always reading one physical block ahead of the program and 
writing one physical block behind the program. When a fatal error 
occurs in a program and the program is outputting to the user's 
console, the output will be interrupted by the monitor error mes- 
sage. Because of the overlapped output and computation, the loca- 
tion of the error in the program usually has no relation to where 
the output was interrupted. 


All I/O is done by physical blocks which are a fixed size for each 
device. The interesting sizes are 128 words for magnetic tapes, 
drums, and disks, and 16 words for the user's console. There is 
no relation between blocks (physical) and records (logical). 


RANDOM ACCESS 1/0 


Random access is implemented for directory devices by means of the 
two calls: 


RREAD(unit,record,size) - random access read 
RWRI (unit,record,size) - random access write 


These calls position read and write pointers respectively within a 
binary file of constant record length "SIZE". After moving the 
pointers, the user may do any I/0 operations which apply to binary 
files. He can read or write a number of consecutive records start- 
ing at the pointers without the overhead of resetting the pointers 
for each record. 


The sophisticated user can use these calls to move the pointers in 
files of non-constant record length. 


APPENDIX 5 


BASIC DIFFERENCES BETWEEN FORTRAN II 
AND ALC FORTRAN IV 


Variable Type 


Variables may be declared by type through the use of the DOUBLE PRECISION, COMPLEX, INTEGER, 
LOGICAL, and REAL type specifications. 'mplicit typing may be accomplished through the use of the 
IMPLICIT specification statement. 


Mixed Mode 

Mixed mode expressions are permitted except for the combination of the double precision and complex 
quantities. | 

Function Naming 


The initial letter of functions is used to type the values of functions. Thus, the LOG, LOGI@, and 
FIX functions have been changed to ALOG, ALOGIG, and IFIX, ete. The terminal F in function 


names is no longer meaningful, and function names may have from one to six characters. 


Arithmetic Function Statement Dummy Arguments 


In FORTRAN IV if a variable appears both as a dummy argument in an arithmetic statement function 


and as an ordinary variable in the same program, its type is the same in both contexts. 


Hardware Tests 
All hardware tests and settings such as IF ACCUMULATOR OVERFLOW and SENSE LIGHT i have been 
changed to subroutine calls such as CALL OVERFL(i) and CALL SLITE(i). 


Input/Output 


The following input/output statements have been changed: 


FORTRAN II FORTRAN IV 
READ TAPE u, list READ — (u)list 
READ INPUT TAPE iu, f, list READ  (u, f)list 
WRITE TAPE u, list WRITE (u,)list 
WRITE OUTPUT TAPE uf, list WRITE = (u, f)list 


a ee 


Y 7. COMMON and EQUIVALENCE 


In FORTRAN IV, EQUIVALENCE does not affect the ordering within common blocks. EQUIVALENCE 
may only have the effect of lengthening a common block. COMMON statements may contain dimen- 


sion information. 


8. EXTERNAL 


Arguments of subprograms which are external subprograms are declared as such through the use of the 


EXTERNAL statement. 
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APPENDIX 6 


ALC FORTRAN IV COMPILER DIAGNOSTICS 


Diagnostic Cause 


ALLOCATION legal dimension specification, COMMON or 
| EQUIVALENCE statement with a dummy 
argument subprogram name or an inconsistency 
in storage assignments in a COMMON or 
EQUIVALENCE statement. 


ALLOCATION ERRORS Inconsistent allocation of variables in COM- 
MON and/or EQUIVALENCE statements. 


CONST OVFLO Too many significant digits in the formation of 
a constant. 

CONTINUATION CARDS More than 19 continuation cards. 

DATA CNT Incorrect number of constants supplied for a 


DATA statement variable list. 


ID CONFLICT Use of scalar where array is required, CALL to 
a scalar or array, subprogram with a CALL to 
itself, attempt to dimension an array more than 


once in one DIMENSION statement. 


ILLEGAL CONVERSION IMPLIED Attempt to mix double precision and complex 7 


data in the same expression. 
ILLEGAL DO CLOSE Illegal statement terminating a DO loop. 


ILLEGAL IF ARG Logical IF or DO statement adjacent to a logical 


IF or illegal expression within a logical IF. 


LABEL IWlegal statement label field. 


16. 


18. 


19. 


21. 


Dia tic 
MULTIPLY DEFINED LBLS 


NO STATEMENT 


NOT ENOUGH SUBSCRIPTS 


NUMBER 
OPEN DO LOOPS 


ORDER — 


SYNTAX 


TYPE 
TOO MANY SUBSCRIPTS 


UNDEFINED LBLS 


UNTERMINATED HOLLERITH STRING 


ey ar 


Cause 


- Twoor more identical statement labels. 


A statement label appears with an empty state- 


ment field. 
An array variable appears with too few subscripts. 


Hlegal statement label as in 
GO TO 999999 


Terminating statement for DO loop(s) omitted or 


improper nesting of DO loops. 


Statements out of order such as COMMON, 
DIMENSION, EQUIVALENCE, TYPE, DATA, 
NAMELIST, SUBROUTINE, FUNCTION, or 
BLOCK DATA appearing after an executable 


statement. 


Compiler cannot recognize statement as a 


properly constructed FORTRAN IV statement. 
Use of noninteger subscript or DO parameter. 
An array variable appears with too many subscripts. 


Statements referenced by GO TO, IF, READ, 
WRITE, etc. have been omitted. 


A missing single quote or fewer than n characters 


following an "nH" specification. 


10 


900 


20 


50 


40 


60 


70 


970 


APPENDIX 7 


SAMPLE PROGRAM 


THIS PROGRAM DOES NOT DO ANYTHING USEFUL 
EXCEPT TO DEMONSTRATE USE OF MANY OF THE 
I/O FEATURES OF THE ALC FORTRAN IV SYSTEM. 


Statement 


DIMENSION IRAY(100) 


CALL IFILE(9,'DATA','SEQ") 


CALL SETSEQ(9) 


READ(9,900)N,A 
FORMAT (1, F) 


IF (N-100) 20,30,30 
DO 40 I=0,24 

DO 50 K=1,100 

TRAY (K) =100*1I+K 
WRITE(3)IRAY 

REWIND 3 

SKIPRECORD 3 
IF(EOFC(3))70,60,60 
CALL RREAD(3,12,100) 
READ (3) IRAY 

TYPE 970,I1RAY(1) 
FORMAT (1X,16) 


CALL RWRI(3,26,100) 


Comment 


Open file "DATA.SEQ" for input 
on Unit 9. 


Set flag to read sequence numbers 
on "DATA.SEQ" as ASCII data. 


Read the sequence number into "N" 
and a floating point number into 
wat 

The horizontal tab in the sequence 
number acts as the illegal char- 
acter for the free field format. 
Exit on sequence number 100. 


Write 25 records of length 100. 


Put consecutive integers in IRAY. 


Write in "FILE3.DAT". 
Close file. 


Skip to end of file. 
Read pointer at record 12. 
Read record 12. 


Type first integer on user console. 


Set write pointer at record 26. 


e 


id 


Do 80 I1=1,3 Write record 12 on records 26, 
27, and 28. 


80 WRITE (3) IRAY 


30 END Enough. 


